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BALLISTIC INTERCEPT MISSIONS TO COMET ENCKE 


FINAL REPORT 
OF THE 

COMET ENCKE BALLISTIC MISSION ENGINEERING PANEL 


1.0 INTRODUCTION 

A ballistic intercept of Comet Encke presents a new and unique challenge in the 
exploration of the solar system. The planetary orbits (with the exception of 
Mercury and Pluto) are all within a few degrees of the ecliptic plane and have 
nearly circular orbits. Comet Encke 's orbit on the other hand, is inclined ap- 
proximately 12° to the ecliptic plane and has an eccentricity of 0. 85. Comet 
activity is a function of the distance between the comet and the sun. These fac- 
tors present an intriguing challenge to determine the optimum ballistic intercept 
of a spacecraft with the comet. If optimum intercept is defined as that mission 
which yields the highest degree of scientific value for the resources expended, 
then the following factors must be considered in the synthesis of a mission. 

Energy requirements to effect an intercept are minimal if the intercept occurs 
approximately 30 days prior to perihelion. However, encounter velocity is high- 
est for this trajectory. Launch energy requirements increase rapidly as the 
intercept occurs closer to perihelion. This increase in energy dictates that the 
launch propulsion requirements increase, or the spacecraft weight and hence 
the science payload decrease. However, an intercept near perihelion results in 
the minimum relative velocity between spacecraft and the comet. As the inter- 
cept point is varied in heliocentric distance, the spacecraft targeting error also 
varies. The combination of encounter conditions, targeting error, and comet 
activity dictate how close the aim point should be to the cometary nucleus. Dis- 
tance to the nucleus at closest approach, encounter velocity, and cometary acti- 
vity at intercept are major constraints on the science value of the mission. The 
relative aspect angle of the spacecraft trajectory varies with heliocentric dis- 
tance. At perihelion the spacecraft will encounter the comet at approximately 
90° phase angle. The phase angle decreases for earlier encounters and ap- 
proaches 0° for an intercept 17 days prior to perihelion. Spacecraft engineer- 
ing requirements and restraints are also a function of the intercept point. An 
intercept near perihelion presents the most severe thermal control requirements 
to the spacecraft design. Also communications range at intercept is a function 
of where the intercept occurs. These factors necessarily influence the space- 
craft complexity, and hence, cost. 

The problem of defining the optimum intercept condition is therefore a difficult 
one in that all of the parameters must be assessed from the standpoint of maxi- 
mizing the science value per dollar. This report documents the analysis 
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conducted by the Encke Ballistic Mission Engineering Panel to arrive at a logical 
basis for determining the most profitable ballistic mission to Comet Encke. 


2.0 BALLISTIC MISSIONS TO ENCKE: MISSION ANALYSIS 

Three different intercept missions of Comet Encke are described in this section. 
Case 1 is for nominal encounter at perihelion which occurs about a day after the 
descending node at a flyby speed of about 8km/s, case 2 is for encounter at 16 
days before perihelion (P-16) at a flyby speed of 18km/s, and case 3 is for en- 
counter at P-30 days at a flyby speed of 27km/s. One advantage of case 3 over 
cases 1 and 2 is that the launch energy requirement is significantly less. On the 
other hand, the flyby speed is smallest for the perihelion intercept and the comet 
is most active there. Finally, the targeting errors are different for each mis- 
sion. Launch and trajectory characteristics are discussed in Appendix I. 

2.1 Mission Descriptions 

The short period Comet Encke during its 1980 apparition offers the fortuitous 
possibility of a short flight time mission with a relatively small flyby speed. 

This happens because the Earth crosses the ascending node of the orbit of Encke 
at such a time that near 180° , a Hohmann type transfer trajectory will yield a 
flyby near the perihelion. The time interval between Earth's crossing the node 
on about August 28, 1980 and the perihelion of Encke at 0. 34 AU on December 7, 
1980 is 101 days. 

Figure 1 is an ecliptic plane projection of the orbit of Encke with the position of 
the Earth indicated at the times of the perihelia in August 1977, December 1980 
and March 1984. A bipolar plot is also shown to illustrate the relative positions 
of the Earth and the comet (Fig. 2). (A bipolar plot is one in which the sun-earth 
line remains fixed at a given apparition and thus the motion of the comet with re- 
spect to the earth is easily seen.) The relative positions of Sun, Earth and Comet 
can be displayed for any apparition by locating the Earth at the moment of peri- 
helion with respect to the vernal equinox. Positions of the Earth at the perihelia 
of 1977, 1980, and 1984 are shown. (The vernal equinox direction indicated is 
valid only at the time of perihelion. ) Table I contains the orbital elements for 
the comet in 1980 (Marsden and Sekanina, (Astron, J. March 1974)). 

Figure 3 shows the relative motion of the spacecraft with respect to the comet 
near encounter on each of the three mission trajectories. Each curve is drawn 
projected in the plane determined by Sun-comet line at encounter and the rela- 
tive velocity vector at encounter. The spacings of tick marks at one day intervals 
along each curve indicate the difference in flyby speeds for the three cases. The 
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Figure 1. Comet Encke Orbit Plot (Ecliptic Plane Projection) 


-100 days 



BELOW 

Figure 2. Comet Encke Bipolar Plot With Earth for 1977, 1980 
and 1984 Apparitions (Ecliptic Plane Projection) 
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Figure 3. Comet Encke Encounter Geometry (1980) 


Table I 

Orbital Elements for P /Encke 


q = 

0.339941 AU 

T p = Dec. 6.5470 1980 

e = 

0.846470 

Period = 3. 295 yr. = 1203. 4 days 

i “ 

11.9460 

Aphelion Distance = 4.0884 AU 

a » 

334.1978 

Speed at q = 69.416km/s 

GO = 

185.9792 

a = 2. 214167 AU 


estimated nuclear magnitude is indicated at several points. It is based on the 
standard inverse square law (for objects shining by reflected light), an absolute 
magnitude of 16.0 and a phase increase of 0.03 mag/deg. (See Section 3. 1.) 

Finally Table n contains mission and trajectory data for each mission including 
the range of values corresponding to a 10 day launch period and a 4 or 5 day arri- 
val period (except for case 1 where a single arrival date is chosen in order to 
achieve the lowest possible flyby speed, see Appendix I). 
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Table II 


Encke 1980 Mission Data 


Case 1: 

Perihelion Encounter 


Launch Vehicle Titan- 3 E/Centaur/TE 364-3 


Range 

Nominal 

Launch Date 

August 20-29 

August 25 

Arrival Date 

December 8 

December 8 (P+2) 

Launch C 3 (km 2 /s 2 ) 

90-104 

94 

Payload System Weight (kg) 

1059-823 

977 

Sun-Comet Distance (at 
Encounter, AU) 

0.34 

0.34 

Earth-Comet Distance (at 
Encounter, AU) 

1.08 

1.08 

Phase Angle (Degrees) 

83-90 

87 

Thermal Load (at Encounter, 
Solar Constants) 

8.7 

8.7 

Flyby Speed (km/sec) 

7. 1-8.2 

7.3 

Case 2 : 

0. 53 AU Encounter 


Launch Vehicles 

(i) Atlas/Centaur/TE 364-4 

(ii) Titan-3E/Centaur 

Launch Date 

August 24 -Sept, 3 

August 28 

Arrival Date 

Nov. 29 - Nov. 24 

November 22 (P-16) 

Launch C 3 (km 2 /s 2 ) 

66-80 

72 

Payload System Weight (kg) (i) 

363-282 

327 

(ii> 

1550-1150 

1400 

Sun-Comet Distance (at 
Encounter, AU) 

0.42-0.56 

0.53 

Earth-Comet Distance (at 
Encounter, AU) 

0.53-0. 72 

0.58 

Phase Angle (Degrees) 

0.5-8 

0.5 
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Table II (Continued) 


Case 2: (Continued) 



Range 

Nominal 

Thermal Load (Solar Constants) 



At Encounter 

3. 2-5. 7 

3.6 

At Perihelion of S/C Trajectory 

4. 3-6. 2 

5.0 

Flyby Speed (km/sec) 

15-19 

18 

Case 3; 

0 . 80 AU Encounter 


Launch Vehicle Atlas/Centaur/TE 364-4 

Launch Date 

July 31-Aug. 10 

August 8 

Arrival Date 

Nov. 4-Nov. 8 

November 6 (P-30) 

Launch C 3 (km 2 /s 2 ) 

41-46 

43 

Payload System Weight (kg) 

580-534 

569 

Sun-Comet Distance (at 
Encounter, AU) 

0.76-0.82 

0.80 

Earth- Comet Distance (at 
Encounter, AU) 

0.30-0.37 

0.34 

Phase Angle (Degrees) 

12-14 

12 

Thermal Loads (Solar Con- 
stants) At Encounter 

1.4-1. 8 

1.6 

At Perihelion of S/C Trajectory 

1.4-2. 4 

1.9 

Flyby Speed (km/sec) 

25-27 

27 


2.2 Navigation 

2.2.1 A-Priori Position Errors for the 1980 Apparition of Comet Encke 

The position errors derived in this section are those which result from uncer- 
tainties in the position of Comet Encke using groundbased observations only. 

For the present covariance error analysis, the 5 returns to perihelion (1967- 
1980) are represented by forty actual observations from August 2, 1967 through 
October 24, 1973 and by 28 additional, postulated observations from October 24, 
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1973 through November 16, 1980. One observation was processed at each of the 
1978 and 1979 opposition dates and the 1980 recovery of the comet was assumed 
to occur on July 9. The postulated observation schedule was determined after 
considering the relative -sun- earth-comet geometries, the available hours of 
dark observing time as well as the apparent nuclear and total magnitudes for 
various dates. 

The present error analysis of comet Encke assumes a l-o measurement error 
of 3 arc seconds for both the right ascension and declination. This value is pri- 
marily due to the deviations of the comet's center of mass from the observed 
center of light. The 3 arc second value is consistent with the mean residuals ob- 
tained for various orbit determinations for past apparitions of comet Encke. 

This value is somewhat higher than for most other short period comets due to 
comet Encke 's relatively high nuclear activity. 

The error analysis was initialized in 1967 by a state vector and appropriate 
values for the nongravitational parameters. The initial 8x8 covariance matrix 
was essentially infinite. Each set of observations was batch processed and 
the updated covariance was propagated forward in time via the state transition 
matrix. All planetary perturbations were taken into account. The time history 
of the comet's error ellipsoid is presented in Table HI. The first two columns 
of this table represent the calendar date and the number of days to perihelion. 

The next twelve columns represent the la position errors (km) for the radial 
Sun-comet direction (r), the direction normal to the comet's orbital plane (n), 
and the transverse direction defined by the cross product of the first two unit 
vectors (T = n x r) . The columns headed by A , R and 6 represent the Earth- 
comet distance in AU, the Sun-comet distance in AU and the Sun-Earth-Comet 
angle in degrees. The a-priori errors represent the forward propagation of the 
covariance matrix obtained by processing all observations from 1967-1979. The 
three columns listed under cases 1, 2 and 3 reflect the effect of each 1980 ob- 
servation on the comet's error ellipsoid. Horizontal lines are drawn under the 
last ground based observation in each of the three cases. For example, for 
case 1, the final ground based observation on November 16 reduces thea r , o R 
and a T components to 416, 249 and 359 km. In the absence of further observa- 
tions, a T grows to 897 km and ff r shrinks to 143 km at the December 8th inter- 
cept. This behavior is due primarily to the dynamic evolution of the error el- 
lipsoid in the absence of observations. From the table, it is apparent that, for 
case 1, the cometary position errors (o T , o n , a r ) at intercept are 143,249 and 
897km. For case 2 the respective values are 422,240 and 450km, and for case 
3, the final position errors at intercept are 422, 295 and 282 km. 

2.2.2 On-Board Optical Navigation 

The use of on-board optical measurements of Encke r s nucleus is being con- 
sidered to reduce the a-priori ephemeris errors (Sec. 2.2. 1), and consequently 


7 



Table IE 


A-Priori Position Errors for Comet Encke (1980) 








1980 Observations Processed* 

* 





Date 1980 

T 

P 

A 

-Priori* 

Case 1 Intercept 
Dec. 8 

Case 2 Intercept 
Nov. 21 

Case 3 Intercept 
Nov. 6 

<AU) 

(DEG) 



a r 


0 J 

ff r 

°n 


a r 

a n 

°T 



0 J 

A 

R 

e 

July 9 

-150 

4168 

2130 

3239 

3352 

1926 

2471 

3352 

1926 

2471 

3352 

1926 

2471 

2.43 

2.33 

72 

19 

-140 

4338 

2084 

3275 

2917 

1739 

2012 

2917 

1737 

2012 

2917 

1737 

2012 

2.21 

2.23 

78 

29 

-130 

4521 

2036 

3327 

2572 

1567 

1683 

2572 

1567 

1683 

2572 

1567 

1683 

1.98 

2.13 

84 

Aug. 8 

-120 

4718 

1985 

3399 

2271 

1406 

1426 

2271 

1406 

1426 

2271 

1406 

1426 

1.76 

2.03 

90 

18 

-no 

4933 

1936 

3494 

1992 

1249 

1213 

1992 

1249 

1213 

1992 

1249 

1213 

1.53 

1.92 

96 

28 

-100 

5169 

1894 

3617 

1644 

1146 

1026 

1644 

1146 

1026 

1644 

1146 

1026 

1.31 

1.81 

101 

Sept. 7 

-90 

5429 

1876 

3769 

1469 

945 

836 

1469 

945 

836 

1469 

945 

836 

1.10 

1.69 

107 

17 

-80 

5717 

1921 

2952 

1217 

799 

710 

1217 

799 

710 

1217 

799 

710 

0.89 

1.56 

111 

27 

-70 

6040 

2117 

4146 

968 

658 

564 

968 

658 

564 

968 

658 

564 

0.69 

1.43 

114 

Oct. 7 

-60 

6403 

2604 

4301 

724 

524 

427 

724 

524 

427 

724 

524 

427 

0.51 

1.28 

112 

17 

-50 

6811 

3480 

4347 

504 

400 

313 

504 

400 

313 

504 

400 

313 

0.36 

1.13 

103 

27 

-40 

7258 

4612 

4407 

387 

308 

264 

387 

308 

264 

387 

308 

264 

0.28 

0.97 

77 

Nov. 6 

-30 

7699 

5595 

5229 

391 

269 

273 

391 

269 

273 

422 

295 

282 

0. 32 

0.80 

45 

16 

-20 

7893 

5683 

8023 

416 

249 

359 

1 416 

249 

359 




0.47 

0.62 

29 



Table III (Continued) 


Date 1980 

T P 

A 

-Priori* 

1980 Observations Processed** 

(AU) 

(DEG) 

Case 1 Intercept 
Dec. 8 

Case 2 Intercept 
Nov. 21 

Case 3 Intercept 
Nov. 6 


a n 

Oj 

a r 

c n 

°T 

CT r 

ff n 

o T 

a r 

°n 

°T 

A 

R 

9 

Nov, 21 

-15 

7617 

5181 

10263 

422 

240 

450 

422 

240 

450 




0.58 

0. 52 

25 

26 

-10 

6628 

4481 

12910 

401 

234 

579 







0. 70 

0.44 

23 

Dec. 6 

0 

399 

3477 

16833 

171 

243 

874 







1.00 

0.23 

20 

8 

+2 

1871 

3411 

16611 

143 

249 

897 







1.06 

0.34 

19 


*A priori, one-sigma errors (Km) in the radial, normal and transverse directions. Last observation 
processes was during 1979 opposition. 

** Evolution of one-sigma errors (Km) if one ground based observation is processed at 10 day intervals begin- 
ning July 9. Measurement noise = 3 arc seconds. 

For cases 1, 2 and 3, the tabular entries just above the horizontal line represent the position errors after 
the last ground based observation. 




to enable improved targeting accuracy. The key to such measurements is locat- 
ing the nucleus relative to an inertial frame, i. e. , a star background. The re- 
quirements which an on-board system must satisfy are: 

a. detection of the nucleus, 

b. detection of sufficient comparison stars, 

c. determination of their positions with sufficient cross track accuracy to 
improve the epheme rides. 

On-board navigation (OBN) should not be confused with the imaging experiment, 
which also features optical detection of the nucleus, OBN requires an attitude 
stabilized camera whereas the imaging experiment could be accomplished with 
either spin-scan or attitude -stabilized cameras. A de-spun platform {such as 
on OSO) would be required for OBN on a spin-stabilized spacecraft. In either 
ease, spacecraft pointing errors could impose special limitations on the maximum 
exposure time. Since all OBN must be done from an attitude stabilized platform, 
we have used the parameters of vidicon cameras already flown on Mariner space- 
craft to assess the possibility of OBN. Other optical sensors may be possible 
{Sec .6.0). 

The threshold magnitude at which Encke's nucleus can be detected by state-of- 
the-art imaging systems has been a somewhat controversial question since it 
requires a detailed appraisal of the technical level of performance of flight- 
proven cameras. A detailed study has been conducted at JPL and shows that a 
10th magnitude detection threshold may be feasible using Mariner-type vidicons. 
The actual thresholds were arrived at in the following way: 

The Mariner telescope (500 mm f. 1.) was assumed to have enhanced UV trans- 
mission due to a change in coatings and filters. The MVM electronic noise 
characteristics were used. The comet nucleus was assumed to be detected when 
the signal to RMS noise was 10. The maximum expo sure times used were 6 seconds 
for the 500 mm telescope (Mariner 9) and 18 seconds for the 1500 mm focal 
length telescope (Mariner 10). Exposure time as a function of magnitude is shown 
in Figure 4. This requires a drift rate of 5 m rad/sec. for 1 pixel smear. Re- 
sults for two albedos were considered (20% and 5%) and a 4 Km nucleus was as- 
sumed in Table IV. These values generally reflect 10th magnitude detection with 
500 mm telescope at pr 2 =0.8 and about 9th magnitude with pr 2 = 0. 2. Use of 
the 1500 mm telescope causes detection to be delayed by 1-2 magnitudes. First 
detection times for a 2 Km nucleus are given in Appendix IV. The time of detec- 
tion before intercept is strongly dependent on the assumed threshold magnitude 
for first detection. 
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Figure 4. Exposure Time vs. Magnitude 
Table TV 


Time of Detection of Encke's Nucleus in Days Before Encounter 



pr 2 

500 mm Telescope 

1500 mm Telescope 



(a) (b) 

(a) (b) 

Mission 1 

0.8 

25 ±5 

12 ±2 

Perihelion Encounter 

0.2 

17 ±4 

7 ±2 

Mission 2 

0.8 

12 ±2 

6.5 ±2 

*= : 0.54AU Encounter 

0.2 

8.5 ±1.5 

4 ±1 

Mission 3 

0.8 

7.5 ±1.5 

4 ±1 

^0,8 AU Encounter 

0.2 

4.5 ±1.5 

2 ±1 


(a) 10th magnitude detection threshold. 

(b) 9th magnitude detection threshold. 
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For the purpose of determining the improvement in a-priori ephemeris errors 
a 9th magnitude detection of the nucleus was assumed (Sec. 2. 2. 3). These re- 
sults are all based on analysis and modelling since demonstration of such de- 
tection has not been possible. Thus, we can say the probability of detecting 
Encke's nucleus as a function of nuclear magnitude is less than unity. Such 
probability must be factored into the prediction of error estimates for the mis- 
sion's navigation. When this is done it is found that it alters significantly the 
notion of one to three sigma error ellipses, tending to lower the probability of 
successful targeting to very small circles around the nucleus. This will be con- 
sidered below in the discussion of targeting errors. 

Star background plots for the three missions being considered are shown in Fig- 
ures 5 through 7. Mission 2 has an abundance of adequate stars; mission 1 ap- 
pears to have an adequate number of stars and mission 3 appears to have insuf- 
ficient stars. The geometry will change slightly during the approach. Stars up 
through 9th magnitude were considered. 

Accuracy turns out not to be an important parameter (over a limited range) in 
that the key requirement is detection. As long as several on-board measure- 
ments can be obtained, improvement in position determination will result. We 
have assumed 30 m rad as the one-sigma measurement accuracy consistent with 
Mariner 9 and Mariner 10 experience. If it is twice as bad, then acceptable per- 
formance - as reported in the next section - will still result. 



Figure 5. Background Stars at Perihelion 
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Figure 7. Background Stars at P-30 
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2.2.3 Targeting Errors 

The targeting errors described in this section are derived from the a-priori 
cometary position errors (Sec. 2.2. 1) and the reduction in these errors which 
results from on-board navigation (Sec . 2. 2. 2) and earth-based video tracking of 
the spacecraft. The approach geometries relative to the sun-comet line for each 
of the three model intercepts are shown in Figure 3. The number of position 
measurements obtained with the on-board imaging system is determined by tak- 
ing measurements at 6 hour intervals, starting when the nucleus reaches 9th 
magnitude as seen by the spacecraft, and stopping when the distance to the nu- 
cleus reaches 2 x 10 6 Ion or at one day prior to encounter, whichever occurs 
sooner. The 9th magnitude limit was discussed in the previous section and it 
is assumed that continuous undisturbed measurements are desirable within 2 x 
10 6 Km of the nucleus to ensure detection of a bow shock, and that system engi- 
neering response times require completion of navigation activities at least one 
day prior to intercept. 

Targeting errors are represented by 3a ellipses in the "impact plane, 11 which is 
normal to the spacecraft velocity vector relative to the comet at encounter. (The 
closest approach relative position vector always lies in this plane.) Figure 8 
shows the "T-R” coordinate system used. "T" is defined by the intersection of 
the impact plane with the ecliptic parallel which passes through the comet at en- 
counter. 6 and the relative velocity vector, V, complete a right-hand system. 
The expected 3a error ellipses for the three encounter dates under consideration 
are illustrated in Figure 9 with the numerical statistics given in Table V. Both 
comet and spacecraft errors are reflected in these ellipses, so the comet posi- 
tion can be considered a random variable about the center with the spacecraft a 
fixed point chosen anywhere in the plane. To minimize miss distance while 
avoiding passing through an exclusion zone about the nucleus an "aim point" 
should be selected co-linear with the minor axis of the ellipse. 

Although a significant improvement of targeting errors would be possible should 
on-hoard na vigation be employed, it should be remembered that the feasibility 
of on-board navigation for comet-intercept missions has not yet been demon- 
strated. For a conservative appraisal of the potential scientific return from the 
1980 Encke flyby, the targeting errors without on-board navigation should be 
used (see Sec. 4.0) . 

2. 3 Targeting Strategy 

The targeting strategy has been developed assuming: (1) that the spacecraft must 
pass on the sunward side of the nucleus in order to ensure favorable conditions 
for the i maging experiment; (2) that there may be some exclusion zone (of radius 
R EZ ) within which the particle impact hazard to the spacecraft is unacceptably 
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V- SPACECRAFT VELOCITY RELATIVE TO COMET 
Z ■ NORTH ECLIPTIC POLE 

A /N A 

T s V « Z ' REFERENCE VECTOR IN IMPACT PLANE 

AAA 

R ■ V X T IN IMPACT PLANE 

= ORIENTATION OF MAJOR AXIS OF ERROR ELLIPSE 
IN IMPACT PLANE 

Figure 8. Impact Plane Coordinate System 




WITHOUT ON-BOARD NAVIGATION WITH ON-BOARD NAVIGATION 

( NOTE SCALE CHANGE ) 

Figure 9. Targeting Error Ellipses (3a) in the Impact Plane Showing 

Three Encounter Dates 
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Table V 


Targeting Errors 


Encounter 

P + 2 Days 

P - 16 Days 

P - 30 Days 

On-Board Navigation 

Without 

With 

Without 

With 

Without 

With 

Error Ellipse (3a) 







Semi-Major Axis (Km) 

1141 

279 

1414 

334 

1003 

208 

Semi-Minor Axis (Km) 

389 

63 

509 

155 

537 

112 

Orientation- 4 > (Degrees) 

73 

60 

18 

123 

63 

159 

On-Board Navigation Assumptions 






Measurement Accuracy (1-a) 

6 arc sec 




Measurement Schedule 


One every 6 hours 



First Measurement 


Point where Encke’s Nucleus ~ 9th Mag. 

Last Measurement 


S/C-Comet Distance ~2 

x 10 6 Km 



high; and (3) that the median miss distance for each heliocentric intercept point 
will depend on whether or not on-board navigation can be used to reduce the 
ephemeris errors. 

The first assumption is satisfied a-priori at the P-30 and P-15 intercepts since 
the spacecraft approaches the comet essentially along the sun-comet line (Fig. 3) 
It is satisfied at perihelion (where the spacecraft approaches at 60-90° to the 
sun-comet line) if the spacecraft is targeted 3o off the nucleus in the sunward 
direction. The second assumption requires that the spacecraft be targeted at 
R EZ = 3c in all three intercepts if Rez =£ 0. R EZ can only be determined by a 
detailed hazard analysis which is beyond the scope of this study; for comparison 
purposes we have taken ad hoc values: R EZ = 400 Km and R EZ = 0. The third 
condition together with the expected scientific return in each case (see Section 
4.0) allows an evaluation of the need for on-board navigation. 

The three targeting geometries are shown in Figure 10. The probability of fly- 
ing by the nucleus at a distance between R and R + dR is given by P(R)dR and 
has been calculated for each of the twelve cases using Monte-Carlo techniques 
(Section 2.2.3). The median miss distance is listed in Table VI for each case. 
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CASE A: TARGETING STRATEGY ASSUMING AN EXCLUSION 


ZONE OF RADIUS R £z . Rf 0.34;053 ; 0- 8 AU 




CASE B: TARGETING STRATEGY WHEN APPROACHING FROM 
SUN ALONG THE SUN - COMET LINE. R EZ ■ O 

R e « 0.8 ; 0.53 AU 



CASE C : TARGETING STRATEGY WHEN APPROACHING AT 

RIGHT ANGLES TO THE SUN-COMET LINE. R EZ “0 


R e = 0.34 AU 



Figure 10 
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Table VI 


Median Miss Distance vs. Intercept Conditions 


R(A. U. ) 

On-Board 

Median Miss Distance (Km) 

Navigation 

R 

EZ 0 

^EZ 

•400 Km 

0.34 

Yes 

95 ^ 

| 

► Case A* 

471 ' 



No 

510 J 


865 


0.53 

Yes 

98 ' 


566 

► Case C 


No 

403 

► Case B 

1017 


0.80 

Yes 

62 


515 



No 

302 j 


996 . 



^Targeting strategy, see Figure 10. 


3.0 COMET ACTIVITY MODEL 

A model of the physical characteristics and activity of Comet Encke is essential 
to enable an estimate of the scientific return for each of the proposed intercept 
missions. Fortunately, Comet Encke has been studied fairly intensively over 
the past 150 years, consequently a relatively larger body of data exists for it 
than for many other short -period comets. Some of these data, such as the trans- 
verse nongravitational acceleration - A 2 , are relatively straightforward and 
error free while others may be influenced by substantial experimental errors, 
as is the case with brightness measurements at low elevations and at small ang- 
ular distances from the sun. We have tried to incorporate the latest knowledge 
of Comet Encke in formulating a baseline model which includes the basic ele- 
ments necessary to estimate the science return for the different missions. 

The major scientific questions to be investigated on a flyby mission deal with the 
nature of the nucleus and the nuclear region (icy halo 7 ) , the identity and abun- 
dances of the parent molecules, the nature and size (and spatial) distribution of 
the dust, the nature of the solar wind interaction with the comet, and the mech- 
anisms leading to destruction and ionization of the cometary gas. The current 
state of our technology determines the kinds of on-board instruments that may 


18 




be used to study a particular problem. Thus the nucleus may be studied by re- 
mote optical scanning (e. g. , by a TV-system) ; the gaseous components by direct 
sampling {ion and neutral mass spectrometers) and indirect sampling (floures- 
cence radiation in characteristic band systems or atomic lines) . The dust com- 
ponent may be studied directly using impact ionization detectors and micro- 
meteoroid impact counters and indirectly by observations of Miescattered 
sunlight. Magnetic and electric field detectors may be used to determine the 
ambient fields directly, and electron energy analyzers will yield data on the 
abundance and energy distribution of free electrons. These data will be neces- 
sary for a definitive understanding of the physics of Comet Encke. 

The following discussion of the physical activity of Comet Encke is organized so 
as to correspond to the proposed studies of nucleus, gas, and dust. 

3. 1 Nucleus 

It is generally agreed that the observed non-gravitational effect on the motion of 
comets is persuasive evidence for the existence of a solid icy nucleus. The non- 
gravitational effect (N.G.E.) is expressed as a vector acceleration (A) with com- 
ponents resolved along the sun-comet line (Aj ) , perpendicular to the sun-comet 
line in the orbit plane (A 2 ) and mutually perpendicular to the first two directions 
A 3 ) (see Fig. 11). The N.G. E. arises when gas evolves asymmetrically from 
the nucleus. If the nucleus is not spinning only A j is non- zero; a spinning nucleus 
can lead to non-zero instantaneous values for A 1 , A 2 , A 3 and non-zero time aver- 
aged values for A l and A 2 . A is more easily and accurately determined than 
Aj since it is observed as a change in the time of perihelion on successive 
returns. 

The values of A 2 for Comet Encke have been derived from orbital data of the 
past 150 years by Marsden and Sekanina and are shown in Figure 12. A 2 in- 
creased from 1800 to 1825 then decreased by more than a factor of ten to a very 
small value by 1971. (Aj has also been reduced from the orbital data (Yeomans, 
private communication) but the resulting accuracy is so poor that the formal 
values are positive as often as negative, which is physically unrealistic.) Sekanina 
has suggested that the decrease in A 2 can be explained by a decrease in the gas 
production rate and that this is consistent with the derived secular decrease on 
the brightness of the comet over this time. However, A 2 may decrease in mag- 
nitude because of secular decreases in the gas_production rate, or precession of 
the nuclear spin axis toward the orbit plane. A 2 is now known to have changed 
sign for two short period comets indicating that their spin axes passed through 
their respective orbital planes (Fig. 12) . 

The photometric brightness of Comet Encke does appear to show some secular 
decrease, although the amount of this decrease is a topic of some dispute (Kresak, 
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A, ALMOST ALWAYS POSITIVE (ENCKE SLIGHTLY NEGATIVE) 

A 2 POSITIVE OR NEGATIVE, MAY SHOW SECULAR DECREASE, 

MAY CHANGE SIGN (ENCKE’S A 2 DECREASING) 

A|/A 2 TYPICALLY ~ 10 BUT << 10 FOR ENCKE 

Figure 11. Nuclear Rotation and Non-Gravitational Acceleration 



Figure 12. Secular Changes in Transverse Non-Gravitational 
Acceleration of Several Short Period Comets 
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Vanysek, Meisel, Sekanina, Bortle). A quantitative determination of the de- 
crease in brightness is made difficult because of greater dust loads in the air 
in modern times, extensive background sky brightness at low elevations, depend- 
ence of magnitude estimates on telescope aperture, the great variability in view- 
ing geometry for successive apparitions, increasing light pollution, and so forth. 
Beyer, in this century, and Holetschek, in the last, are responsible for the two 
sets of observations which are most internally consistent and which span several 
apparitions. Beyer's magnitude estimates for successive apparitions have been 
reduced to a common magnitude law by Sekanina by introducing residual magni- 
tudes (mi). Sekanina derived a secular decrease of 0. 1 magnitudes per revolu- 
tion from Beyer's data. We have performed a similar renormalization and find 
somewhat smaller values for the residual magnitudes and thus for the secular 
decrease (< 0.0 5 magnitudes per revolution). Beyer's data (renormalized by 
Mumma) and data taken during the 1974 apparition are shown in Figure 13. 

Bortle (private communication), using the same data but a slightly different anal- 
ysis, concludes that Encke shows no secular decrease. It should be noted that 
our residual magnitudes are essentially constant for 1937 and 1947, and again 
nearly constant for 1951 through 1974. Thus, except for a step function between 
1947 and 1951 the brightness would seem to be nearly constant. (Mumma be- 
lieves this step function is the result of increased sky brightness over Hamburg 
after World War n, tending to lower the apparent cometary brightness.) How- 
ever, for the purposes of this report, we adopt Sekanina's value for the secular 
decrease (0. 1 magnitudes per revolution) since this establishes a lower limit 
to the expected total brightness in 1980. 

The brightness clearly tends toward an inverse square dependence (n=2) on 
heliocentric distance (R) at small R, which is consistent with the assumption 
that the available solar radiation limits the gas production rate. This is made 
reasonable by the following argument. The intensity of light radiated by a 
particular molecular species is given by 

Ii ~ 8* Nj (1) 

where I. is the total number of photons radiated per second into 4 ir steradians by 
the species in question, g| is the solar g-factor for the molecular band system 
or atomic line in question, and Nj is the total number of the species in the coma 
in the lower electronic state. Also, 

N i ~ Qi T i (2) 

where Qj is the production rate and t { is the mean lifetime, or 

l i ~ 8A T i ~ (§iVi R ’ 2 Q i r2 (3) 
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Figure 13. Total Brightness vs. Heliocentric Distance 


or 


Ij-teiViQi (4) 

If Qj is limited by the available solar radiation then 

Ii ~ (gi^iQi)i R " 2 (5) 
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i. e. , the limiting brightness law is as the inverse square of R. For the past 12 
returns to perihelion, Encke's brightness has regularly approached an inverse 
square law (R<0. 8 AU), in agreement with Equation (5). 

The total magnitude (T mag ) is given by 

T mag = 8.9 + 5 log A + 5 log R (6) 

from R — 0 . 8 AU to R - 0. 34 AU (perihelion) . Equation (6) is the relation obtained 
if all subsequent observations are reduced to the 1937 apparition. The corres- 
ponding brightness for the 1980 return will be smaller by the secular decrease 
between 1937 and 1980 or by about 1.4 magnitudes if Sekanina's estimate of 0. 1 
magnitudes per revolution is accepted, i. e. , 

T mag = 10.3 + 5 log A + 5 log R 0.34 AU < R < 0.8 AU (7) 

We also adopt the following law for the nuclear magnitude (N mag ) based on 
Roemer's data (see Appendix H), 

N mag = 1 6.0 + 5 log A + 5 log R + 0.03(3 (8) 

where (3 is the phase angle, A the comet-observer distance in AU, and R is the 
sun-comet distance in AU. The nuclear magnitude is assumed to remain in- 
variant through 1980. In contrast to Encke's regular behavior prior to perihelion, 
the comet sometimes shows flaring and irregular behavior several weeks after 
perihelion (Dossin, Roemer) . Furthermore, the brightness after perihelion is 
generally about a factor of three (or one astronomical magnitude) weaker than 
at corresponding distances before perihelion (Sekanina) . This behavior was ob- 
served to hold in 1974 (Bennett) when Encke's brightness at P + 30 days was found 
to be about one magnitude smaller than that at P - 30 (Fig. 13) . Thus, while 
Encke appears to fade somewhat well after perihelion, the pre-perihelion data 
show that Encke's brightness is developing as expected at least to within one day 
of perihelion. 

Delsemme has shown that reasonable estimates of the nuclear radius and Bond 
albedo can be obtained for comets Tago-Sato-Kosaka (1969 IX) and Bennett (1970 
II) by comparing the gas production rate at small heliocentric distances and the 
brightness at large heliocentric distances where reflected sunlight dominates the 
brightness. Unfortunately, this procedure does not work well for Comet Encke 
and Delsemme interprets this to mean that the nucleus of Encke is only partially 
ice-covered or that the vaporization of water ice does not control the production 
rate of hydrogen. However, a lower limit to the radius (r min ) may be established 
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by assuming that all of the solar radiation at 0.7 AU is effective in vaporizing 
water ice, i.e., 


F S R " 2 » - Qh 2 o l < 9 > 

where F Q = solar flux at 1 AU (0.03 cal cm -2 sec -1 ), R = heliocentric distance 
(AU), Qh 9 o = production rate of H 2 0 molecules/sec, and L is the heat of subli- 
mation of water ice near 200 K (11,500 cal/mole or 1.9 x 10' 20 cal/molecule) . 
Then, 


r min 


0.17Km 


( 10 ) 


Other limits on the size of the nucleus can be established from the measured 
brightness at distances large enough that sun light reflected from the nucleus is 
the dominant source of brightness. In this case, 

log [pD 2 0(0)] = 0.4 (M q - M c ) + 2 log RA (11a) 


where p is the geometric albedo, D is the nuclear radius, 0(0) is the phase law, 
0 is the phase angle, M 0 = solar magnitude, and M c the cometary magnitude. 
This can be written 


N mag = M c = M o +5k)g rA + ai0 (lib) 

where a l - 2. 5 0' 1 log [0(0)1 . Roemer’s data give pD 2 = 0. 7 (taken at aphelion 
in 1972) and pD 2 -0.13 (at r = 2. 9 AU in 1973) (Appendix II) . Observations in 1963 
give pD 2 —0.22 and aj = 0.013. A geometric albedo of ~0.1 would indicate a 
radius of ~2. 7 Km using the 1972 data, although a wide range of effective radii 
can be derived depending on the assumed albedo. An upper limit of ~6 Km is set 
by assuming a geometric albedo of 0. 02 which is less than the lowest known 
asteroidal albedo (0.027). 

The import size parameters for the purpose of this study are the minimum 
radius (0. 17 Km) and the nominal design radius (chosen to be 1 Km) which to- 
gether allow us to model the science return of the imaging experiment (Sec. 4.1.2). 

3, 2 Gas 

As discussed earlier, insolation appears to limit the production of gas for helio- 
centric distances less than about 0. 8 AU. We therefore take as our baseline 
model for parent molecule production a rate which increases as R -2 . We nor- 
malize the H 2 0 production rate at 0. 71 AU to one-half the hydrogen atom produc- 
tion rate observed by Bertaux et al. , in 1971 (the adopted secular decrease would 
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further reduce this by about 25% in 1980) . Thus, 


^-,(0.7AU) = 2.5 X 10 26 steradiaiT 1 sec' 1 (12a) 

4tt 


and 


^ H2 ° (R) = 1.2 X 10 26 R' 2 steradian" 1 sec _1 ,RinAU (12b) 

4n 

We also assume uniform gas production per cm 2 of area from a spherical nucleus 
(point source model) . This leads to lower gas densities in the sunward direction 
than if we had assumed emission only from the sunward hemisphere, but higher 
densities than would be present within an icy halo. A simple flow model then 
leads to gas densities in the coma given by 


n(D, R) 


n 0 (l km, 1 AU) v Q 
"r 2 v t D 2 



(13) 


where V T is the terminal flow velocity and t 1 is the lifetime of the parent mol- 
ecule against photodissociation at R = 1 AU (Fig. 14). It did not seem worthwhile 
to use the frictional flow model of Mendis, Holzer, and Axford because the gas 
densities and corresponding collision frequencies are low and because a spheri- 
cal model is only an approximation anyhow. The initial thermal velocity (v D ) at 
the surface of the nucleus changes slowly with heliocentric distance. The rela- 
tive production rates of all other parent molecules are taken to be constant rela- 
tive to water, such as would be the case for a clathrate icy nucleus. For con- 
venience trace molecules are each assumed to be 1% as abundant as H 2 0. 


The photodissociation lifetime of H,0 at 1 AU is taken to be ~7 x 10 4 sec and is 
assumed to vary as R 2 . All other species are assumed to have the lifetimes 
given by Potter and del Ducca, which are subject to some uncertainty but are the 
only ones available at the present time. 


3. 3 Dust 


The dust distribution has been modelled after Finson and Probstein's results for 
Arend-Roland, and the brightness of Mie-scattered sunlight from dust particles 
in the field-of-view of the imaging experiment has been determined (Appendix 
III). The dustscattered background light (watts/cm 2 srpm) is 10 2 to 10 3 times 
smaller than the light reflected from the nucleus showing that the imaging experi- 
ment should have no trouble 'seeing' the nucleus through the intervening dust. 

The large difference between the nuclear brightness and dust brightness allows 
for a sizeable error in the dust abundance. Even a 10-fold increase in dust 
would not prejudice imaging of the nucleus. 
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Figure 14. H 2 0 Density vs. Distance from the Nucleus 


A detailed dust-impact hazard analysis is beyond the scope of this study. 


4.0 EXPECTED SCIENCE RETURN 

4. 1 Experiments Which Are Strongly Dependent on the Mission Mode 

4.1.1 Neutral Mass Spectrometer 

The instantaneous counting rate of a molecular species seen by a neutral mass 
spectrometer at encounter is given by 

dN _ 

— = EG n(D, R) (14) 

at 

Where E is the detection efficiency (taken to be 1 count/sec per 1000 molecules/ 
cm 3 ), G is the duty cycle, and n is the local number density at a distance D from 
the nucleus (Eq. 13) . The number of counts detected for a particular ion per 


26 




unit kilometer of flight path is given by 

dN _ _j dN 
dL “ V dt 

or 


(15) 


dN 

dL 


GE Q(R) 
4tt(L 2 +D2 i )V t 


exp[-(L 2 +D 2 n ) ,/ VV T r 1 R 2 ] 


= C1(L 2 + D 2 ) exp[-C2(L 2 + D 2 ,)*] 


(16) 


dN/dL has been calculated for each mission concept for 16 different sets of data, 
representing the combinations of the values for (Q, r, D M ) given in Table VII. 

The entries labeled 1% Qh 9 o apply either for (G = l, 1% H 2 0 abundance) or for 
(G = 0. 01, H 2 0 abundance) . The results for dN/dL are shown as the solid lines 
in Figures 16 through 31. 


Table VII 


Baseline Data for Neutral Mass Spectrometer Study 



— 1 — — ■ '■ . 

Mission 1 

Mission 2 

Mission 3 


0.34AU 

0.53 AU 

0.80 AU 

C1: (Qh 2 o) 

3.2E8 

6.1E7 

1.8E7 

(1% Qh 2 o) 

3.2E6 

6.1E5 

1.8E5 

C2: (rj = 7E4 sec) 

6.9E5 

2.8E-5 

1.25E-5 

(Tj =3. 5E3) 

1.38E-3 

5.6E-4 

2.5E-4 

D m :R EZ =0, w-OBN 

95 

98 

62 

wo-OBN 

510 

403 

302 

R ez =400, w-OBN 

471 

566 

515 

wo-OBN 

865 

1017 

996 
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The number of measurements of a species on a given flyby has been estimated by 
integrating from Lj to L 2 along the flight path (beginning at 32, 000 Km) until the 
error in the accumulated counts (N -1/2 ) is equaled or exceeded by the error in 
cometocentric distance (AD/D) for that measurement (Fig. 15). An 'X' is then 
plotted at the center of the interval (Lj ,L 2 ) and the process is begun again at 
L 2 . The distance D at which the combined measurement error falls below 10% 
is marked by a vertical line on each curve. The number of measurements ob- 
tained having errors less than 10% on the flight path in to D M is shown at the 
upper left of each curve, and the % error of the measurement at D = D M is also 
shown. In some cases, the % error of the first measurement is shown at the 
extreme lower right. The science return can thus be compared for various mis- 
sions by inspection. 



Figure 15. Flyby Geometry at Encounter 
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A note of caution should be introduced. If a sizeable amount of H 2 0 is locked in 
icy grains in the inner coma, the H,0 densities will begin to fall progressively 
further belowthe values given by the point source model as D decreases from the 
value (sayD H ) corresponding to the outer limits of the icy halo (D H ~ 1000 Km?). 
Any departure of the gas densities from a D“ 2 distribution would only be recog- 
nized on those intercept missions for which the accuracy had already reached 
10% or better at distances greater than D H . 

The first case considers measurements of H 2 0 with 100% duty cycle and shows 
that all intercept missions return adequate measurements of H 2 0 (Figs. 16-19). 
Even for the worst case considered (no ori-board navigation, and a 400 Km ex- 
clusion zone), H 2 0 measurements reached 10% accuracy at ~4,000 Km from the 
nucleus at 0. 8 AU and ~74 measurements can be expected. The minimum dis- 
tance of approach is ~1000 Km in all three missions, but at 0, 34 AU measure- 
ments begin at 20,000 Km and~260 measurements are obtained. An icy halo may 
never be entered in these cases, however. 

The case of a trace constituent with the lifetime of H 2 0 (~20R 2 hours) but only 
1% as abundant is shown in Figures 20-23. This can be viewed as the science 
return for H 2 0 if a 1% duty cycle mass spectrometer is used. Then the 10% 
error point is reached before entering the icy halo (1000 Km) only on the peri- 
helion intercept. On-board navigation would definitely be required to return 
adequate data. 

The next case considers a hypothetical major eonsituent but with a one-hour life- 
time at 1 AU (Figs. 24-27). The production rate is assumed equal to the H 2 0 
production rate. The 10% error point is reached at nearly the same distance 
(D ~ 2000 Km) on all three missions, in spite of the short lifetime. The science 
return is adequate for all cases except that for R EZ = 400, and without on-board 
navigation (Fig. 27) for which the icy halo is barely penetrated; however, the 
number of measurements even for this case is adequate outside D = 100 Km. 

The final case considers a trace molecule 1% as abundant as H 2 0 but with a one- 
hour lifetime at 1 AU (Figs. 28-31) . Only a perihelion intercept with R EZ = 0 and 
on-board navigation (Fig. 28) has a hope of adequate data return and then only if 
the radius of an icy halo is small (say ~100 Km) or if the fraction of trace vola- 
tiles locked an ice grains is small. 

4.1.2 Imaging 

A detailed study of the imaging experiment has been carried out. We have as- 
sumed a vidicon camera of the Mariner-1971 B-type on an attitude stabilized 
platform (either 3 -axis or de-spin). The brightnesses of the coma and nucleus 
are as described in Section 3.1 and in Appendix IV. The miss distances of 
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Section 2. 3 have been used in modelling the imaging results for each of the 12 
missions treated (with/without optical navigation; with/without an exclusion zone 
for R = 0. 34, 0.53, 0. 80 AU) . The results of interest are: total number of pic- 
tures returned, spatial extent per picture element (pixel) for the best frame, 
and the number of days before encounter for earliest nucleus detection (this im- 
pacts on the feasibility of terminal navigation) . 

The results are summarized in Tables VIII and DC, for further details see 
Appendix IV. 

4. 1. 3 Optical Particle Detector; Impact Detector and Analyser 

It is assumed that an optical particle detector canbe built, patterned after Sisyphus 
but without the obvious problems encountered by Sisyphus. It is further as- 
sumed that with proper background rejection filtering that the instrument will 
have a "range to particle radius ratio" of at least 10 5 at 0. 8 AU for geometric 
albedos of 0. 1. For calculation purposes the nominal symmetric dust model of 
Taylor, et al, , is assumed. No account is taken of any change of emission with 
heliocentric distance. Inside 1000 km (cometocentric distance) additional counts 
can be expected from icy grains. Because of the great uncertainty in size and 
number of these grains, they are not included in the comparison, but they ob- 
viously will increase the count on 100 and 500 km flybys (also see Appendix IH). 


Table Vffl 

Comet-Sun Encounter Range Effects on Imaging 


Parameter 

Case I 

Case n 

Case m 

Encounter-Sun Distance (AU) 

0.34 

0.53 

0.80 

Flyby Velocity (km/s) 

8.6 

18.0 

26.8 

Encounter Duration (sec) 
2 km > 5 pixel s 

3000 

1480 

545 

Smear Rate at Encounter 
(pixels /sec) 

287 

600 

895 

Earliest Nucleus* 
Resolution (days) 

E - 3 

E - 1.6 

E - 0.8 

With 10% coma filter 

E - 7 

E - 2.5 

E - 1.6 


♦Assumes: M = 16 + 51ogR + Slog A + 0.03g 

Nucleus Diameter = 2 km 
p r 2 = 0.24 
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Table IX 

Comet-Sun Encounter Range Effects on Imaging 


Parameter 

Case I 

Case n 

Case III 

50 K bits/sec 

120 K bits/sec 

120 K bits/sec 

No. of Pictures* w/o Recorder 




Large Frames (84 sec/fr) 

33 (4**) 

18 

7 

Large Frames (10 sec/fr) 

33 (4**) 

39 

14 

Small Frames (15 sec/fr) 

94 (8**) 

39 

36 

No. of Pictures* w Recorder 




Large Frames 

36 

18 

7 

Large Frames 

240 

148 

55 

Small Frames 

200 

148 

55 


♦Assumes: Large frame = 4.5x 10 6 bits; 

Small frame = 1.6 x 10 6 bits; 

Despun Platform 
♦♦Without despun platform; 

Bit rate is 4 K bits/sec 

There are only two large differences in the three mission classes. The particles 
emitted are brighter nearer the Sun, the "range to radius ratio" increasing lin- 
early with solar distance, and the count rate at any given distance from the comet 
is obviously a function of the spacecraft velocity. The tables and graphs for this 
instrument therefore indicate the total counts as a function of size and the count 
rate at closest approach to the cometary nucleus. 

When the spacecraft velocity is much higher than the particle velocities it can 
be assumed simply that the spacecraft is flying through a static distribution of 
debris. The particle concentration C and the total particle fluence during the 
flythru F are given by 


E 

C = — and F - 

4 7T R 2 ^ 


E 

4V 


where E = particle emission rate from nucleus (#/s) 
R = distance from nucleus (m) 
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v = particle velocity (m s' 1 ) 

R 0 = distance of closest approach to nucleus (m) 

The optical particle detector looks out in a fixed direction with a half cone angle 
a to a range r determined by the particle size, The maximum number of parti- 
cles per second N and the total number of counts T are given by 

N = C R V r 2 tan a 

, 0 

and 

T ~ F r 2 tan a 

where V = velocity of the spacecraft (ms -1 ). The equation for T is quite inac- 
curate for larger particles, since they can be seen at ranges which are large 
relative to the distance to the nucleus , distances at which the particle concen- 
trations are very different from those in the immediate vicinity of the space- 
craft, The actual values would be more difficult to calculate and would be a 
function of the instrument look angles. 

It is assumed that the impact detector has an area of 10 " 2 m 2 . Everything scales 
linearly with detector area, if some other value in fact pertains. The total num- 
ber of counts t will just be the fluence on that area A. The maximum count rate 
n is just the concentration at closest approach, times the detector area, times 
the spacecraft velocity. 

t = FA n - CVA 


The data on dust taken from Taylor, et al. , is as follows: 


Size Range 
(m) 

Emission Rate (E) 
(s' 1 ) 

Velocity* (r>) 
(m s' 1 ) 

Calculation Size (r ) 
(m) P 

3 x 10' 7 

- 10' 6 

2.5 x 10 11 

300 

6.5 x IO' 7 

IO" 6 

- 10 - 5 

9 c; v ia12 

u ( v A J.V 

1 CO 
IdU 

5 x IO" 6 

10' 5 

- 10' 4 

2.5 x IO 10 

50 

5 x 10' 5 

io- 4 

- 10' 3 

2.5 x 10 7 

15 

5 x 10' 4 

10' 3 

- 10' 2 

2.5 x 10 3 

3 

5 x 10' 3 

10" 2 

- 10" 1 

2.5 x 10" 1 

0.1 

5 x 10' 2 

10' 1 

- 2. 5 x 10' 1 2.5 x 10~ 2 

0.01 

1.75 x 10' 1 


•These approximate velocities of the particles were calculated for 1 AU heliocentric distance and would be 
somewhat larger nearer the Sun. 
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0 . 8 AU 26 Km/ sec 100 km approach 


Size (m) 


C (m 3 ) 

F <m" 2 ) 

N (s' 1 ) 

T 

3 x 10" 7 - 

IO" 6 

6,63 x IO -3 

2.08 x IO 3 

1.3 x IO" 1 

1. 5 x 10° 

IO -6 - 

10' 5 

1. 33 x KT 1 

4.17 x 10 4 

1. 5 x 10 2 

1. 8 x IO 3 

10" 5 - 

10-4 

3.98 x 10' 3 

1.25 x IO 3 

4.6 x 10 2 

5. 5 x IO 3 

M 

o 

1 

10-3 

1.33 x 10“ 5 

4.17 x 10° 

1.5 x 10 2 

1. 8 x IO 3 

10' 3 - 

IO' 2 

6.63 x 10“ 9 

2.08 x IO" 3 

7. 6 x 10° 

9.2 x IO 1 

IO' 2 - 

10- 1 

1.99 x IO' 11 

6.25 x IO' 6 

2.3 x 10° 

2.8 x IO 1 

IO' 1 - 

2,5 x 10-1 

1.99 x IO' 11 

6.25 x IO' 6 

2.8 x IO 1 

3.4 x IO 2 



0.8AU 26 Km/sec 

500 km approach 


Size (m) 


C (m -3 ) 

F (m- 2 ) 

N (s- 1 ) 

T 

3 x 10~ 7 - 

l 0 -6 

2. 65 x IO -4 

4.17 x 10 2 

5.1 x IO" 3 

3.1 x IO" 1 

10 ~ 6 - 

IO' 5 

5.31 x IO' 3 

8.33 x IO 3 

6. lx 10° 

3. 7 x 10 2 

IO' 5 - 

10- 4 

1.59 x 10- 4 

2.50 x 10 2 

1.8 x IO 1 

1.1 x IO 3 

io- 4 - 

IO' 3 

5.31 x IO" 7 

8.33 x IO' 1 

6.1 x 10° 

3.7 x IO 2 

IO" 3 - 

10- 2 

2,65 x IO' 10 

4. 17 x IO" 4 

3.0 x IO' 1 

1. 8 x IO 1 

10' 2 - 

IO' 1 

7.96 x IO" 13 

1.25 x IO' 6 

9. lx IO -2 

5. 5 x 10° 

10' 1 - 

2, 5 x 10 -1 

7.96 x IO' 13 

1.25 x IO" 6 

1.1 x 10° 

6.7 x IO 1 


0. 8 AU 26 Km/sec 1500 km approach 


Size (m) 


C (m- 3 ) 

F (m -2 ) 

N (s-*) 

T 

3 x 10~ 7 - 

IO" 6 

2.95 x IO- 5 

1. 39 x 10 2 

5.70 x IO" 4 

1.03 x 10-1 

IO" 6 - 

10- 5 

5. 89 x IO -4 

2.78 x IO 3 

6.74 x IO” 1 

1.22 x IO 2 

10 " 5 - 

10 - 4 

1.77 x 10~ 5 

8.33 x IO 1 

2.02 x 10° 

3.67 x 10 2 

IO -4 - 

10-3 

5.89 x IO" 8 

2.78 x IO" 1 

6.74 x IO' 1 

1.22 x 10 2 

IO' 3 - 

10- 2 

2.95 x IO" 11 

1.39 x 10~ 4 

3.37 x IO' 2 

6.11 x 10° 

w 

O 

l 

IO' 1 

8.84 x IO" 14 

4.17 x IO' 7 

1.01 x 10“ 2 

1.83 x 10° 

H 1 

O 

1 

2.5 x 10' 1 

8. 84 x IO -14 

4. 17 x IO' 7 

1. 24 x IO' 1 

2.25 x IO 1 
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0. 53 AU 18 Km/sec 100 km approach 


Size (m) 

C <m -3 ) 

F (m -2 ) 

N (s' 1 ) 

T 

3 x 10~ 7 - 10' 6 

6.63 x 10' 3 

2.08 x 10 3 

2.02 x 10" 1 

3.53 x 10° 

10 -6 - 10" 5 

1.33 x 10 _1 

4.17 x 10 4 

2.39 x 10 2 

4.18 x 10 3 

10" s - 10 -4 

3.98 x 10" 3 

1.25 x 10 3 

7. 18 x 10 2 

1.25 x 10 4 

10" 4 - 10" 3 

1. 33 x 10 -5 

4.17 x 10 2 

2.39 x 10 2 

4.18 x 10 3 

lO" 3 - 10" 2 

6.63 x 10 “ 9 

2.08 x 10 -3 

1.20 x 10 1 

2.09 x 10 2 

10' 2 - lO* 1 

1.99 x 10" 11 

6.25 x 10“ 6 

3, 59 x 10° 

6.26 x 10 1 

10 _1 - 2. 5 x 10" 

1 1. 99 x 10” 1 1 

6.25 x 10" 6 

4.40 x 10 1 

7.67 x 10 2 


0.53AU 18 Km/sec 

500 km approach 


Size <m) 

C (m -3 ) 

F (m~ 2 ) 

N (s" 1 ) 

T 

3 x 10" 7 - 10" 6 

2.65 x 10~ 4 

4.17 x 10 2 

8.09 x 10" 3 

7.06 x 10" 

10" 6 - lO" 5 

5.31 x 10" 3 

8.33 x 10 3 

9. 57 x 10° 

8.35 x 10 2 

10 " 5 - lO" 4 

1.59 x 10" 4 

2.50 x 10 2 

2. 87 x 10 1 

2.51 x 10 3 

10 " 4 - 10" 3 

5.31 x 10" 7 

8.33 x 10 _1 

9.57 x 10° 

8.35 x 10 2 

10 " 3 - 10" 2 

2.65 x lO" 10 

4.17 x 10" 4 

4. 79 x 10" 1 

4.18 x 10 1 

10" 2 - 10 -1 

7.96 x 10“ 13 

1.25 x 10" 6 

1.44 x 10" 1 

1. 25 x 10 1 

10" 1 - 2.5 x 10" 

1 7.96 xlO" 1 3 

1. 25 x 10“ 6 

1.76 x 10° 

1.53 x 10 2 


0.53AU 18 Km/sec 1500km approach 


Size (m) 

C (m* 3 ) 

F {m" 2 ) 

N (s _1 ) 

T 

3 x 10" 7 

- 10" 6 

2. 95 x lO -5 

1.39 x 10 2 

8. 99 x 10~ 4 

2.35 x 10” 

10~ 6 

- 10 - 5 

5. 89 x 10~ 4 

2. 78 x 10 3 

1.06 x 10° 

2.78 x 10 2 

10" s 

- 10- 4 

1. 77 x 10 -5 

8.33 x 10 1 

3.19 x 10° 

8.35 x 10 2 

1 

O 

i— \ 

- 10- 3 

5.89 x lO’ 8 

2.78 x 10 _1 

1.06 x 10° 

2.78 x 10 2 

10" 3 

- 10- 2 

2.95 x 10 -11 

1.39 x 10" 4 

5.32 x 10" 2 

1. 39 x 10 1 

10" 2 

- 10" 1 

8.84 x 10“ 14 

4. 17 x 10- 7 

1.60 x 10~ 2 

4. 18 x 10° 

I - 1 

O 

1 

- 2.5 x 10' 1 

8.84 x 10" 14 

4. 17 x 10- 7 

1,95 x 10" 1 

5.12 x 10 1 
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0 . 34 AU 8 Km/sec 100 Ion approach 


Size (m) 

C (m -3 ) 

F (m -2 ) 

N (s _1 ) 

T 

3 x 10“ 7 - 10" 6 

6.63 x 10~ 3 

2.08 x IO 3 

2.18 x IO" 1 

8.58 x 10° 

IO’ 6 - 10“ 5 

1.33 x IO" 1 

4.17 x 10 4 

2.58 x IO 2 

1.01 x 10 4 

IO" 5 - IO" 4 

3.98 x IO" 3 

1.25 x IO 3 

7. 75 x IO 2 

3.04 x 10 4 

IO -4 - 10~ 3 

1.33 x 10 ~ 5 

4.17 x 10° 

2. 58 x 10 2 

1.01 x 10 4 

10 -3 - 10“ 2 

6. 63 x 10 " 9 

2.08 x IO -3 

1.29 x IO 1 

5.07 x 10 2 

10~ 2 - IO -1 

1.99 x 10 " n 

6.25 x IO" 6 

3.88 x 10 2 

1.52 x 10 2 

IO -1 - 2.5 x 10 

_1 1.99 x 10 “ n 

6.25 x IO" 6 

4.75 x IO 1 

1.86 x IO 3 


0 . 34 AU 8 Km/sec 

500 km approach 


Size (m) 

C (m- 3 ) 

F (m -2 ) 

N (s _1 ) 

T 

3 x 10- 7 - 10- 6 

2.65 x IO" 4 

4,17 x IO 2 

8.74 x IO" 3 

1.72 x 10° 

10~ 6 - IO” 5 

5.31 x IO" 3 

8.33 x IO 3 

1.03 x 10* 

2.03 x IO 3 

10" 5 - 10~ 4 

1.59 x IO" 4 

2. 50 x 10 2 

3.10 x 10 l 

6.09 x 10 3 

10" 4 ~ IO" 3 

5.31 x 10“ 7 

8.33 x IO" 1 

1.03 x IO 1 

2.03 x IO 3 

10~ 3 - IO -2 

2.65 x IO" 10 

4. 17 x IO -4 

5. 17 x IO" 1 

1.01 x 10 2 

lO- 2 - 10' 1 

7.96 x 10~ 13 

1.25 x 10~ 6 

1. 55 x IO" 1 

3.04 x IO 1 

10 _1 - 2.5 x 10- 

■* 7.96 x IO" 1 3 

1.25 x 10~ 6 

1. 90 x 10° 

3. 73 x IO 2 


0.34 AU 8 Km/sec 

1500 km approach 


Size (m) 

C (m -3 ) 

F (nr 2 ) 

N (s' 1 ) 

T 

3 x 10‘ 7 - IO" 6 

2.95 x IO" 5 

1. 39 x 10 2 

9. 71 x IO” 4 

5. 72 x IO" 1 

10-6 _ 10‘ 5 

5. 89 x IO" 4 

2.78 x 10 3 

1.15 x 10° 

6.77 x 10 2 

IO -5 - lO" 4 

1.77 x 10“ 5 

8.33 x IO 1 

3.45 x 10° 

2.03 x 10 3 

IO " 4 - io - 3 

5. 89 x 10“ 8 

2.78 x IO’ 1 

1.15 x 10° 

6.77 x 10 2 

IO" 3 - 10“ 2 

2.95 x IO" 11 

1.39 x 10“ 4 

5. 74 x 10" 2 

3.38 x IO 1 

10' 2 - 10 _1 

8.84 x 10~ 14 

4.17 x 10~ 7 

1.72 x IO -2 

1.01 x IO 1 

10' 1 - 2.5 x 10' 

1 8.84 xlO -14 

4.17 x 10“ 7 

2.11 x IO" 1 

1. 24 x 10 2 
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100 km approach 


Size (m) 

26-km/s 

18km/s 

8km/s 

3 x 10" 7 - 10' 6 

1.72 x 10° 

1.19 x 10° 

5.30 x 10" 1 

10 -6 - 10" s 

3.46 x 10 1 

2. 39 x 10 1 

1.06 x 10° 

10 ‘ 5 - 10" 4 

1.03 x 10° 

7.16 x 10 _1 

3.18 x 10 _1 

10 “ 4 - 10" 3 

3. 46 x 10" 3 

2.39 x lO" 3 

1.06 x 10‘ 3 

10 " 3 - 10 -2 

1.72 x 10" 6 

1.19 x 10' 6 

5. 30 x 10" 7 

10~ 2 - 10" 1 

5.17 x 10" 9 

3.58 x 10" 9 

1.59 x 10“ 9 

10 _1 - 2.5 x 10 _I 

5.17 x 10" 9 

3.58 x 10~ 9 

1. 59 x 10" 9 


500 km approach 


Size (m) 

26km/s 

18km/s 

8km/s 

3 x 10“ 7 

- 10- 6 

6. 89 x lO' 2 

4.77 x 10“ 2 

2.12 x lO -2 

10 ~ 6 

- 10“ 5 

1.38 x 10° 

9.56 x 10” 1 

4.25 x lO"' 

10" s 

-10- 4 

4. 13 x 10" 2 

2.86 x lO' 2 

1.27 x 10~ 2 

10 " 4 

O 

1 

w 

1. 38 x 10" 4 

9.56 x 10- 5 

4.25 x 10~ 5 

1— » 
O 
1 

- 10- 2 

6.89 x 10~ 8 

4.77 x 10“ 8 

2. 12 x 10~ 8 

10" 2 

- 10- 1 

2.07 x lO -10 

1. 43 x lO" 10 

6.37 x 10" 11 

O 

l 

- 2.5 x 10" 1 

2.07 x lO" 10 

1.43 x lO -10 

6.37 x 10' 11 


1500 km approach 


Size (m) 

26km/s 

18km/s 

8km/s 

3 x 10~ 7 

- 10" 6 

7.67 x 10" 3 

5.31 x 10“ 3 

2.36 x 10“ 3 

10" 6 

- 10" 5 

1.53 x 10' 1 

1.06 x 10-* 

4.71 x 10 " 2 

10" 5 

- 10“ 4 

4.60 x 10" 3 

3. 19 x 10- 3 

1.42 x 10“ 3 

1 

O 

rH 

1 

w 

O 

i 

1.53 x 10' 5 

1.06 x lO -5 

4.71 x 10“ 6 

10“ 3 

- 10~ 2 

7. 67 x 10“ 9 

5.31 x 10“ 9 

2.36 x 10 “ 9 

lO" 2 

- 10- 1 

2.30 x 10 -11 

1. 59 x 10 -1 1 

7.07 x 10' 12 

10" 1 

- 2.5 x 10' 1 

2.30 x 10" 11 

1.59 x 10" 11 

7.07 x 10" 12 
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RATE (S 


T 



IN EACH FAMILY OF THREE CURVES, 

THE HIGHER COUNT RATES CORRESPOND 
TO THE CLOSER ENCOUNTERS 


— — 0.34 AU, 8 km /sec 
— 0.53 AU, 18 km /sec 
— 0. 8 AU, 26 km / sec 


Figure 32. Optical Particle Detector (I) Total Counts 
During Encounter 



0.34 AU, 8 km /sec 
0.53 AU, 18 km /sec 
0.8 AU, 26 km /sec 


Figure 33. Optical Particle Detector (N) Maximum Count 
Rate During Encounter 




Figure 34. Impact Detector and 
Analyzer (t) Total Impacts 
During Encounter 



26 km /sec, 0.8 AU 
18 km / sec, 0.53 AU 
8 km /sec, 0.34 AU 


Figure 35. .Impact Detector and Analyzer (n) Maximum 
Count Rate During Encounter 
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4.2 Experiments Which Are Independent of the Mission Mode 

In addition to the experiments already discussed for which science return is 
strongly dependent on the heliocentric distance at intercept, there are some ex- 
periments which are expected to yield satisfactory results regardless of inter- 
cept point. This is not to say that the amount of data returned is independent of 
perihelion distance at intercept, rather that scientific return will be adequate 
for all cases. 

The ion mass spectrometer experiment is expected to yield satisfactory results 
for all intercepts because of the high detection efficiency for ions (100% - as 
compared with 0. 1% for neutral molecules in the neutral mass spectrometer ex- 
periment) . However, because the energy distribution of cometary ions is ex- 
pected to change drastically as the spacecraft moves from the transition zone 
to the inner coma, the focussing properties of the ion mass spectrometer must 
be carefully designed to ensure that the measurements can be correctly related 
to the in-situ ion densities and their velocity distributions (see also Section 6. 3). 
If the design problems can be solved satisfactorily, the information return is 
expected to be adequate for each of the three intercepts discussed here. 

The magnetic field experiments (A.C. and D. C.) and electric field experiments 
(A.C. and D. C.) are expected to return adequate data at all three intercept dis- 
tances. Existing instruments have been used successfully on both three-axis- 
stabilized and spinning spacecraft to measure the fields in the un-shocked and 
shocked solar wind. Since these fields increase with decreasing heliocentric 
distance, the measurements should become progressively easier (e.g. , requir- 
ing progressively shorter antenna length in the case of the E -field experiment) 
as R decreases from 0.8 to 0.34AU. 


5.0 MAJOR ISSUES AS A FUNCTION OF INTERCEPT 

This study has resulted in identifying the following tradeoffs which have major 
impact on the cost and/or scientific return of a ballistic mission to Comet Encke 
in 1980. 


1. A Titan-3E/ Centaur launch vehicle is required for intercept missions 
at 0. 34 and 0. 53 AU whereas an intercept at 0. 8 AU (for a Pioneer S/C 
or one of similar weight) can be accomplished with the much-cheaper 
Atlas/Centaur/TE364-4. A Pioneer can actually be launched to inter- 
cept at 0. 53 AU using an Atlas/Centaur/TE364-4 but the S/C cannot 

handle the thermal loads imposed at this heliocentric distance (see 
(4) on the following page. 
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2. The S/C velocity relative to the comet is lowest (~ 8 Km/sec) for an 
0.34 AU intercept and rapidly increases with heliocentric distance 
(~15 Km/sec at 0.53 and ~26 Km/sec at 0. 8 AU). As we shall see, the 
science return is a strong function of the flyby speed, with higher sci- 
ence return at lower velocities. 

3. The median miss distances (assuming an exclusion zone 400 Kin in rad- 
ius) are approximately the same for all three missions, being ~500 Km 
if on-board navigation is available and ~1000 Km if it is not. If no ex- 
clusion zone is assumed, then the median miss distances are in the 
neighborhood of ~100 Km and ~400 Km, respectively (see Table VI) . 

4. The solar heat load on the S/C is greatest at 0.34 AU presenting an 
engineering problem which will require using the thermal technology 
already developed for the Helios-probe. A Mariner class S/C (e.g. , 
MVM-73) could reasonably be expected to survive and function to 0.4 AU, 
while a Pioneer (e.g. , Pioneer- Venus, if built) could be modified to 
survive to ~0. 7 AU. 

5. The S/C approaches the comet at a small angle to the sun-comet line 

in the intercepts at 0.53 and 0. 8 AU, but at a steep angle (60-90°) to the 
sun-comet line at 0. 34 AU (Fig. 3) . The small angle approach geometry 
has the advantage of traversing the coma and part of the tail with a 
single S/C, but imaging after closest approach is degraded since the 
S/C must look sunward to photograph the nucleus. In addition, separa- 
tion of spatial -temporal variations will be difficult. The high angle 
approach geometry allows one to choose the phase angle at closest ap- 
proach. Also, a 'fish-hook' trajectory is possible so that two bow-shock 
crossings can be made. In this event a second tail probe is desirable 
(launched on the same rocket as the coma probe) and will enable cor- 
related simultaneous measurements on the coma and tail during the 
cometary encounter and in interplanetary space during the cruise mode. 

6. The images of the nucleus (or nuclear region) will have highest resolu- 
tion (~6 m, in 0° phase angle is chosen) and be most numerous (up to 
240) for the 0. 34 AU intercept. The best resolution will be ~55 m (phase 
angle = 90°) at 0. 53 AU and -150 pictures will be returned; at 0. 80 AU 
(90° phase) the best resolution will be —150 mand — 55 pictures will be 
returned. The number of pictures quoted here is twice the number 
which can be taken on the way in to encounter. It may not be possible 

to get all of the images after encounter in the 0. 53 AU and 0. 8 AU mis- 
sions since the camera must slew toward the sun leading to increased 
scattered sunlight. However, it is the close approach imaging which is 
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of greatest value, consequently, loss of some pictures after closest 
approach is not too important. The expected imaging science is 
superior on the 0. 34 AU mission. 

7. An attitude stabilized (3-axis or de-spun) camera platform will be re- 
quired if on-board navigation is deemed desirable. Terminal guidance 
is feasible on all three missions and would be needed to ensure the tar- 
geting accuracy needed to guarantee entry into the icy halo (say R H 
•-'1000 Km) . 

8. Measurements of neutral parent molecules are strongly dependent on 
the heliocentric distance at intercept because of varying cometary activ- 
ity (as R" 2 ) and spacecraft flythrough velocity (as v _1 ). On the other 
hand, the lifetime of molecules in the solar radiation field goes as R 2 . 

A detailed model of the neutral mass spectrometer experiment shows 
that the science return is greatest by far at 0. 34 AU and is roughly in 
the ratio (14/5/1) at (0. 34/0.53/0. 8 AU) respectively. If a mass spec- 
trometer can be built which measures all mass peaks simultaneously 
or is dedicated to a single mass peak (100% duty cycle M. S.), then ade- 
quate data of major constitutents (H 2 0) can be returned on all missions. 
However, adequate data on trace molecules (1% as abundant as H 2 0) 
will be obtained only near perihelion and using on-board navigation. If 
a scanning mass spectrometer (say a duty cycle of 1%) is used, then 
even measurements of H 2 0 are only adequate with on-board navigation 
and at 0. 34 AU. Trace molecules could not be measured with the latter 
instrument in any mission studied. It should be noted that the instru- 
mental sensitivity used in this model could possibly be increased by 
5-10 (A.O. Nier, private communication) with a consequent increased 
science return; and that the model assumes a stabilized platform for 
the neutral mass spectrometer experiment. 

9. Various other desirable experiments have been briefly reviewed and it 
is thought that currently space-proven hardware would return adequate 
data on all missions hut would be maximum at 0.34 AU. 


6.0 RECOMMENDATIONS FOR FURTHER STUDY 

1 . Several experiment problem areas require further study and/or develop- 
ment work which should begin immediately to assure completion/ 
availability in time for the 1980 mission. They are: 

a. Charge-couple-devices and their applicability for use in spin-scan 
imaging systems. 
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b. Development of a 100% duty-cycle cometary neutral mass 
spectrometer. 

c. A detailed study of the ion-mass spectrometer design problem for 
cometary flyby missions. 

d. Long lead-time development of computer software and procedures 
for real-time iteration of the terminal navigation targeting problem. 

e. Phase A studies should be carried out for the 0.34 and 0. 53 AU 
intercept missions. The Phase A study of the 0. 34 AU intercept 
should be carried out jointly with ESRO, if at all possible, since 
Helios thermal technology will be needed. 
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APPENDIX I 


1980 ENCKE FAST-FLYBY 

HELIOCENTRIC LAUNCH AND TRAJECTORY CHARACTERISTICS 



APPENDIX I 


1980 ENCKE FAST-FLYBY 

HELIOCENTRIC LAUNCH AND TRAJECTORY CHARACTERISTICS 


In this appendix, parametric information is presented on the characteristics of 
fast-flyby trajectories to comet Encke during the 1980 apparition. The appropri- 
ate launch period for such missions is near the time when Earth passes through 
the orbit plane of Encke on August 28, 1980. It is fortuitous that Encke's posi- 
tion in orbit at this time allows for a nearly minimum energy transfer to inter- 
cept Encke near perihelion. 

Launch energy contours as a function of launch and arrival date for ballistic 
transfers during the 1980 apparition of Comet Encke are shown in Figure 1-1. 

A range of launch dates near the Encke nodal crossing and a range of arrival 
dates both before and after Encke perihelion are shown. Launch and arrival 
dates are indicated in the Julian calendar for convenience, but conversion to the 
Gregorian calendar is easy by noting that the nodal crossing occurs on August 28, 
1980 (Julian day 2444480) and Encke perihelion occurs at December 6, 1980 
(Julian day 2444580), Various launch energy contours from 65 to 120 km 2 /sec 2 
are shown. Trajectories arriving before perihelion are generally Type I and 
those arriv ing after perihelion are Type H. Also indicated on the right-hand 
edge of the figure is the Encke-Sun distance in AU as a function of arrival 
date. The launch energy capability of a given launch vehicle depends upon the 
characteristics of that vehicle and the mass of the spacecraft being launched; 
however, it can be seen from Figure 1-1 that arrival 15 days prior to perihelion 
requires a launch energy of about 74 km 2 /sec 2 for a 10 -day launch window. In 
general, launch energy increases as the encounter date moves nearer to peri- 
helion. 

Of additional concern is the launch declination of the Earth departure trajec- 
tories. This parameter affects range safety and Earth orbital coast time re- 
quirements. Contours of constant departure declination in the range of launch 
and arrival dates described above are shown in Figure 1-2. Launches from KSC 
are generally limited by a launch azimuth of 180° . Higher azimuths are possible 
by expensive dog-leg maneuvers during ascent to orbit, but are generally not 
done. A launch azimuth of 108° implies a launch declination limit of 33° for 
in-plane launches. Higher declinations are possible by expensive plane-change 
maneuvers,' but are also generally avoided. It can be seen from Figure 1-2 that, 
except for a very small region near perihelion, the launch declinations are always 
below 30° ; and hence, no problem exists for this mission. 
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An important encounter parameter is the relative velocity between the space- 
craft and the comet. Contours of constant arrival speed are shown in Figure 1-3 
for the same range of launch and arrival dates shown in Figure 1-1. The rela- 
tive velocities range from about 20 km/sec for arrivals 20 days prior to perihel- 
ion to about 8km/sec for arrivals near perihelion, Pre -perihelion arrival dates 
have arrival velocities which are relatively constant across a 10 -day launch win- 
dow, while a considerable variation is noted for near -perihelion arrival dates. 
However, careful selection of arrival date/launch date combinations for near- 
perihelion arrival dates can provide a 10 -day launch window with an approach 
velocity of about 8 km/sec. 

An equally important trajectory parameter from the point of view of the science 
value of the mission is the direction relative to the Encke-Sun line from which 
the approach to Encke is made. Contours of constant approach phase angle 
(spaceeraft-Encke-Sun angle) are shown over the launch/arrival date window as 
before in Figure 1-4. It can be seen that arrivals at Encke about 10-20 days 
prior to perhelion result in an approach from nearly the direction of the Sun, 
that is, wdth low' phase angle. The direction of encounter becomes more oblique 
for encounter dates near Encke perihelion. Encounters near perihelion are nearly 
at right angles to the Encke-Sun line; and hence, Encke' s tail. 

The combination of approach velocity and phase angle of the approach dictates the 
time spent by the spacecraft within the comet tail. Figure 1-5 indicates this 
parameter of scientific interest. Contours of constant time spent within 20, 000 km 
of the comet tail (the estimated size of the tail) are shown over the launch/arrival 
date window. The maximum time occurs for encounters about 16 days prior to 
perihelion. This maximum is about 10 hours. The time spent in the tail de- 
creases to about 1 to 1.5 hours as the encounter date approaches Encke 
perihelion. 
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Roemer has measured Encke's nuclear magnitude in 1963 and again in 1972 near 
aphelion. We have derived values of pR 2 0(/3) from her nuclear magnitudes using 
the relationship 

log [pR 2 0(j3)] = 0.4(M @ -M C )+ 21og(rA) (11-1) 

where p = geometric albedo, <p( j3) = phase law, M 0 = solar photographic magni- 
tude (-26,41), and M c = Encke's nuclear magnitude. Roemer's data and our 
calculated results are shown in Table IX— 3 . By inspection of Figure n-1 one sees 
that pR 2 0(o) (i. e. , the nuclear brightness) was greater by about three at aphelion 
in 1972 vs. 1963, Yet we know that the total brightness (at r < 1. 2 AU) shows a 
slow secular decrease at about 0. 05 - 0. 1 mag/Rev, i.e. , the 1974 total bright- 
ness was less than the 1964 brightness (0, 87 -* 0. 75, r < 1. 2 AU) . Dossin ob- 
served that Encke was unusually faint ~30 days after perihelion passage in 1964. 
Thus the observational evidence show that the nucleus was unusually faint before 
perihelion (at large r) , but the total gas production rate was normal until some- 
time after perihelion when it fell off more rapidly than usual. 

The phase law can be derived from the data in Table II-l by the following method. 
We want to cast the nuclear magnitude law in the form 

M c = M 0 + 5 log rA + 0 (II-2) 


From Equation II-l, 

M @ = M 0 + 2.5 log [pR 2 ] , r = A = 1,0 - 0, 0(o) = 1 
log [pR 2 0((3)) = 0.4(M o + 2.5 log (pR 2 ) - M c ) + 2 log rA 


so 


M c - M 0 + 5 log rA - 2.5 log 0(j3) 
M c = M Q + 5 log rA + a x /3 


or 


= ,2 5 log 

0 


(H-3) 
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Figure II-l shows log [pR 2 ^] vs. 0. The value of aj =+0.0130 is determined 
from the 1963 data. Since <fi( o) = 1 by definition, the values pR 2 may be deter- 
mined from the intercepts of the lines. They are pR“ = 0.22 (1963), pR =0.70 
(1972) and pR 2 = 0. 12 (1973). The asteroidal phase law is shown for comparison. 
This is the law recommended by Marsden (1AUC 2547) and is the law adopted for 
the purpose of the main study. 


Table II- 1 

Nuclear Brightness (M 2 ) and pR 2 for Encke's Comet 



r 

A 

0 

m 2 

log [pR 2 <P(0) 

pR 2 0(0) 

9/24/63 Ref. 1 

3.11 

2.11 

2.5 

20.2 

- 0.663 

0.216 

9/25/63 

3.10 

2.11 

2.6 

20.2 

- 0.666 

0.216 

10/12/63 

3.00 

2.07 

8.5 

20.2 

- 0.715 

0.192 

12/16/63 

2.51 

2.46 

22.8 

20.3 

- 0.756 

0.175 

perihelion —5/64 







8/30/64 

1.65 

1.08 

36.4 

19.0 

- 1.316 

0.048 

8/15/72 Ref. 2 

4.09 

3.14 

5.5 

20.5 

- 0.194 

0.64 

Aphelion —9/3/ 64 

4.10 

3.09 

0.8 

N.B.R.* 



9/13/72 

4.10 

3.11 

3.1 

N.B.R.* 



9/21/73 (2586) 

2.89 

1.90 

3.8 

-20.5 

— 0.932 

0.116 


1. E. Roemer and R. E. Lloyd, Astronom. J. 71, 443 (1966). 

2. E. Roemer, IAUC 2435 

3. R. E. McCrosky and C. W. Shao, IAUC 2446 

4. E. Roemer, IAUC 2446 

5. E. Roemer, IAUC 2586 
* Now being reduced 
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Figure H-l. Comet Encke: Nuclear Brightness vs. Phase Angle 


II-3 



LMSC/D4 30482 


APPENDIX III 

OPTICAL IMAGE OF A COMETARY NUCLEUS: 
1980 FLYBY OF COMET ENCKE 


W. C. Wells 
R. S. Benson 
A, D. Anderson 
G. Gal 


August 1974 
Contract NAS5-20566 


Radiation Physics Laboratory 
Palo Alto Research Laboratory 
Lockheed Missiles and Space Company 
Subsidiary of Lockheed Aircraft Corporation 
Palo Alto, California 94304 



APPENDIX in 


OPTICAL IMAGE OF A COMETARY NUCLEUS: 
1980 FLYBY OF COMET ENCKE 


I INTRODUCTION 

This is a study of the feasibility of obtaining optical images of a cometary nucleus 
via a flyby of Comet Encke. This study is based on a physical model of the dust 
cloud surrounding the nucleus. Development of the model is based on available 
physical data and theoretical knowledge of cometary physics. Using this model, 
calculations are made of the absolute surface brightness of the dust in the line 
of sight of the on-board camera and the relative surface brightness of the dust 
compared to the nucleus. The brightness is calculated as a function of heliocen- 
tric distance and for different phase angles (sun-comet-spacecraft angle). 

The study of the feasibility of obtaining optical images of the nucleus of Comet 
Encke via a flyby was divided into two parts. First a physical model of the dust 
cloud surrounding the nucleus was derived and then our Mie scattering code was 
employed to calculate the absolute brightness of the dust in the line of sight of 
the on-board camera. The surface brightness of the comet was compared to this 
background brightness along with contributions from debris and ice particles. 
Conclusions were then drawn in light of these results regarding the optimum im- 
aging system for a flyby mission. 


II BACKGROUND OF FLYBY MISSIONS 

At a symposium on the exploration of space held in Washington, D. C. in April 
1959, Whipple (1959) pointed to the possibility of sending a space probe through 
the neighborhood of a comet. Two years later. Swings (1962) gave an extended 
survey on the scientific objectives and feasibility of such a mission at a sympos- 
ium held in Pasadena, California in August 1961. He summarized the existing 
knowledge about comets and pointed to those problems in the physics of the com- 
ets for which investigations by means of space probes would be of special value. 

There are quite a few problems which would be handled best with instruments 
installed in a cometary probe (Lust, 1969) : 

a. The structure and composition of the nucleus, its surface temperature 
and its color. 
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b. Chemical processes related to the formation of the observed radicals 
which take place in the vicinity of the nucleus where the density is high, 
and the ejection velocity of the different particles. 

c. The composition of the eometary atmosphere, the mechanisms of dis- 
sociation and ionization, the size and nature of the dust grains. 

d. The density gradient and partial densities of the different constituents in 
the coma and in the tail. 

e. The interaction with the solar wind and with interplanetary magnetic 
fields, the structure of shock fronts and related problems of plasma phys- 
ics and magneto-hydrodynamics. 

Among the questions which have to be solved in the preparation of a cometary 
probe is which comet should be selected for a first space mission. Much work 
has been put into this problem by different groups in the U. S. A. and in Europe, 
and all these groups came to realize that the difficulties of such a mission are 
far greater than had first been anticipated . It is evident that a periodic comet 
which has been observed for many apparitions and whose orbit has been calcu- 
lated with some accuracy is an easier target than a new comet which suddenly 
appears without being predicted and whose path has to be calculated from a few 
observations taken in small time intervals. But there are other reasons why a 
periodic comet can be reached more easily by a space vehicle than a new one. 

The orbits of almost all periodic comets have relatively small inclination to the 
ecliptic plane, and a motion in the same sense as Jupiter and all other planets 
(direct orbit). The inclinations of the orbital planes of the new comets are, on 
the other hand, randomly distributed between 0 and 180° , so that on the average 
these comets spend only a small fraction of their time close to the ecliptic plane 
while they are passing through the nodes. This causes a severe restriction for 
the launch window, while retrograde motions (inclination between 90° and 180 ) 
cause very large approach velocities of the spacecraft with respect to the comet. 
Expensive terminal guidance would therefore be necessary to force the trajectory 
of the spacecraft into or near to the orbital plane of the new comet and to dimin- 
ish the relative velocity to a slow flyby or to achieve a rendezvous. 

Because of these difficulties, the new comets have been excluded as a first aim 
for a cometary probe by all groups who have investigated the feasibility of such 
a mission in spite of their high scientific interest. 

Different groups in the U.S.A. as well as in Europe have carried out detailed 
feasibility studies during the last years. In 1961, a space mission to a comet 
was considered for the first time in NASA contracted studies of the Scientific 
and Technical Laboratories Inc. (Corben, 1962) . 
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There exist about 50 periodic comets which have been observed during more than 
one perihelion passage, and an additional 30 comets with periods less than 20 
years observed during one apparition only. About 5-7 of these comets appear on 
the average per year. A great amount of work and extensive calculations have 
been carried out by the different study groups to select among these comets the 
objects being fitted for a first mission. Several selection criteria had to be ap- 
plied with respect to the following points: 

a. The position of the comet at the time of intercept must be known with 
very high accuracy, and it became more and more clear that this condi- 
tion is a major constraint for a cometary mission. It seems desirable 
that the probe should approach within 1000 km. If costly midcourse cor- 
rections are to be avoided, the position of the comet must be predictable 
at the time of launch, that is several months before intercept, with this 
accuracy. 

To secure the necessary accuracy in the calculation of the orbital ele- 
ments, it is necessary to start with a well-known orbit calculated from 
previous apparitions and to apply corrections by means of new observa- 
tions made after recovery. Perturbations by planets, especially by 
Jupiter, can be the reason for considerable changes of the orbit. Also 
secular perturbations caused by a mass loss of the comet and non- 
gravitational forces play a role in this respect. Especially the location 
of the comet in its orbit, that is the daily motion, is affected by such 
perturbations. Since small errors in the daily motion will add up and 
result in large errors in the time of perihelion passage, the accuracy of 
this parameter is very important (Porter et al., 1965). It seemed, 
therefore, necessary to exclude all comets which had been observed dur- 
ing one apparition only, and also some of the remaining comets whose 
predicted orbits were not well enough established. 

b. In order to guarantee observations of the comet’s position necessary for 
orbit corrections, the comet should be recovered at least 2 months be- 
fore launch, and it should be seen against a dark sky for several hours 
per night. This imposes the condition that the comet must be at least of 
20th magnitude— the limit for recovery of an object whose position is 
known— several months before perihelion passage. The average bright- 
ness of periodic comets is so small that not many of them fulfill this 
condition. 

The question of the accuracy of cometary orbits has been investigated by 
the groups in the U. S.A. as well as the European ESRO group. The con- 
clusion is that if every effort is made to recover a comet very early, 
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that is about 20th magnitude, it should be possible to correct the pre- 
dicted orbits to the accuracynecessary fora space mission for a number of 
comets. 

c. The comet should be bright enough at intercept to get photometric data 
anrl good quality spectra from ground-based observatories to supplement 
the spacecraft data. This means that the comet should be of 12th magni- 
tude or brighter at the time of perihelion passage, and it should be seen 
against a dark sky for some hours. Since a great percentage of the peri- 
odic comets never become so bright, and since some of the brighter ones 
are in unfavorable positions on the sky near perihelion (close to the sun), 
this is a further restriction of the number of feasible comets. 

d. The ideal launch velocity should not exceed ~16km/sec. This limit is 
chosen rather arbitrarily, it represents for instance the volocity required 
for a two year flight to Jupiter. It turned out, however, that this require- 
ment does not impose a severe restriction, because almost every comet 
that was brighter than 12th magnitude at intercept fulfilled this condition, 

e. Of two otherwise equal missions, that which leads to a smaller relative 
approach velocity between spacecraft and comet and therefore allows for 
a longer stay of the spacecraft in the vicinity of the comet should have 
priority. 

f. All the points mentioned before refer more or less to the technical side 
of the mission. When choosing a comet, one should of course ask which 
mission promises the best scientific results. Comets are of very differ- 
ent activity according to their n age, " that means according to the number 
of perihelion passages they have already made, because with every ap- 
proach to the sun they lose a considerable amount of volatile material, 
and slowly get exhausted and inactive. 

Summarizing the different points, a feasible comet should fulfill the following 
conditions : 

1. A reliable predicted orbit should be available. 

2. The comet should be of at least 20th magnitude 2 months before launch, 
and at that time be well observable for at least 2 hours per night. 

3. It should be of at least 12th magnitude at intercept, and distinctly visible 
from the earth. 

4. The launch velocity should not exceed ~16 km/sec. 
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5. The relative approach velocity between comet and spacecraft should be 
low. 

6. The comet should be interesting from a scientific point of view. 

A mission to Comet Encke in 1980 meets these requirements and should be 
undertaken. 

HI PHYSICAL PROPERTIES 

With an estimate for the size of the nucleus of Comet Encke, and using the icy- 
conglomerate model, we estimate the distribution of non-volatile dust particles 
near the nucleus. 

\ 

A. Nucleus 

Although Whipple's (1950) model for comet nuclei is not universally accepted 
(e.g. , Lyttleton, 1972), it has been used to successfully explain a number of 
cometary phenomena. Not the least of these is the non-gravitational force re- 
quired to reproduce the orbital motion of Comet Encke and several other comets 
(Marsden and Sekanina, 1974 and earlier papers in this series). The non- 
gravitational acceleration derived by Marsden and Sekanina (1974) for Comet 
Encke shows a regular decrease which is consistent with the existence of a por- 
ous, rocky core, within which "dirty ice" is embedded (Sekanina, 1969b). Since 
the gravitational acceleration is proportional to the fractional rate of change of 
mass (as well as the degree of anisotropy of the ejection), it is necessary to esti- 
mate the mass of the nucleus in order to obtain the rate of mass loss itself. The 
mass (or mean density) and radius are also required in order to calculate the 
terminal velocities of ejected dust particles, since expanding vapor from the 
nucleus must lift the particles against gravitational attraction (Delsemme and 
Miller, 1971). Marsden and Sekanina (1971) estimate that Comet Encke loses 
0.03 percent of its mass during each orbital revolution. A spherical nucleus of 
radius R n = 1.7km and mean density lg cm -3 has a mass ~2 x 10 1 6 g, implying 
AM ~6 x 10 12 g per revolution; this mass loss is comparable to the estimates of 
Sekanina (1969a), which are ~10 13 g. If the geometric albedo of the nucleus is 
a n = 0.1, then a n R n 2 = 0.3km 2 which is to be compared with observed values 
which range from 0.24 (Roemer, 1966) to 5)0.82 km 2 , which we derive from 
Roemer's (1972) observation of the apparent magnitude of Comet Encke at aphel- 
ion (m ~ 20. 5). These lower limits result from our assumption that the comet 
was observed at opposition, thereby giving minimum earth-comet distance and 
phase angle. These results are consistent with the recent treatment of Delsemme 
and Rud (1973). 
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B. Non-Volatile Particles 


1. Size Distribution 

Lacking a determination of particle ejection for Comet Encke itself, we use the 
results of Finson and Probstein (1968) for Comet Arend-Roland and those of 
Sekanina and Miller (1973) for Comet Bennet. In particular, we use the distri- 
bution of particle sizes deduced by Finson and Probstein (1968), which they state 
is well -determined over the range 4 < pd <40 pm g cm“ 3 , where p is the bulk 
density of a particle of diameter d. We use the Sekanina and Miller (1973) dis- 
tribution for the smallest particles (down to pd = 0.9pm g cm -3 ). Specifically, 
denoting by n(a) da the number of particles with radii in the range a to a+da, 


Kjla/ajj) 5 (a - l)/a 0 l<a<6.7pm 

K 2 (a/a 0 )- 3 

n(a) = 

K 3 (a/a o r 4 

K 4 (a/a 0 )" 5 a >44.4 pm 


6.7 < a < 14.4 pm 
14.4 <a < 44.4 pm 


The distribution for larger particles fits smoothly onto the distribution of inter- 
planetary particles as given by Whipple (1968) : n~ a -5 . These are shown in 
Figure III-l. The constant a 0 = 1 pm, and the K's are determined by the re- 
quirement that 


(47rp/3) 



a 3 n(a)da = 1 g 


The rate of ejection of particles of radius a is then given by 


P(a) = n(a) 


dM 

dt 


where dM/dt is the total mass-loss from the nucleus in the form of dust. We 
have 


Kj = 8.916 X 10 14 
K 2 = 1.137 X 10 14 
K 3 = 1.642 X 10 15 
K 4 = 7.298 X 10 16 
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Figure m-1. Particle Size Distribution. 
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The particle mass density is assumed to be p = 0.45 g cm -3 , which is ap- 
proximately the value adopted by Whipple (1967, 1968) for meteroids of mass 
m > 10 " 6 g. Super- Schmidt data for the Taurid meteor shower give a mean 
value of about 0.28 g cm -3 (Vemiani, 1969). A somewhat higher density may 
be appropriate for the smaller particles: Whipple (1968) uses p - lg cm -3 for 
m < 10 -6 g (a < 175 pm), and radio data for sporadic and shower meteors give 
p — 0.8g cm -3 (Verniani, 1973). The particle distribution would not be greatly 
changed if we were to use p - 1 g cm -3 , except that n(a) would peak at a - 0.6 pm 
and extend to about 0.55pm. These small particles would make only a small 
contribution to the total scattered light. 

It is of interest to compare our distribution to one derived by Taylor et al. , 

(1973) in a similar manner. This is done in Figure m-2. The disparities arise 
from the different densities and size cutoffs used. 

2. Production Rate 

As discussed above, non -gravitational acceleration of Comet Encke implies a 
loss of ~10 1 3 g per orbital revolution. It is not known how much of this is in 
the form of dust, but Whipple (1967) suggests that the abundance of Taurid mete- 
ors is evidence that perhaps an old comet such as Encke releases a larger frac- 
tion of material in dust than does a "new" comet (see also Delsemme's 1973 
review). Whipple (1967) goes on to estimate that, averaged over one period, 
Comet Encke contributes as much as 3. 5 x 10 6 g sec -1 of meteoritic material 
to the interplanetary cloud. If true, this would imply a much smaller degree of 
anisotropy of ejection than has been determined for this comet (Marsden and 
Sekanina, 1971, p. 1143). The distribution of particle sizes might differ from 
that for "new” comets, having fewer small particles. 

We conservatively assume that the total mass lost per orbit is 10 1 3 g, of which 
10 percent is dust. This gives a simple average mass -loss rate of dM/dt — 

9.6 x 10 3 g sec -1 as dust. As an approximation to the evaporation rate of water 
snow (Delsemme and Miller, 1971; Marsden. Sekanina and Yeomans, 1973), the 
mass-loss varies as r' 2 for r< r c amd is zero at greater distances. The 
mass-loss rate at any particular heliocentric distance is then easily calculated; 
Table m-1 gives the mass-loss at r = 0. 8 AU for various assumed r c . 

Thus, for this simple approximation, the mass-loss rate is not sensitive to the 
value of r c , and dM/dt — 5 xlO 4 g sec -1 is a typical value for r = 0. 8 AU. 

3. Velocity Distribution 

The evaporation of cometary ices (primarily H 2 0) produces a flow from the 
nucleus which accelerates dust particles to their terminal ejection velocities 
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Figure ni- 2 . Comparison of Distribution of Particle Sizes 
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Table HI-1 


Dust Mass-Loss at 0.8AU 


r c 

(AU) 

dM/dt 

(10 4 g sec -1 ) 

1.00 

5.93 

1.25 

5.47 

1.50 

5.17 

2.00 

4.78 

2.808 

4.42 


within ~10 nuclear radii (Probstein, 1968; Delsemme and Miller, (1971). 
Delsemme and Miller (1971) calculate the largest particle which the expanding 
vapor can lift against gravity; for p - 0 . 45 g cm' 3 , R n = 1. 7 km and p n - 1 g cm 
we have 


a max = 3-14 X 10" 17 Z cm 

where the rate of evaporation, Z (molecule cm -2 sec -1 ), depends on heliocentric 
distance. In Figure in-3 we show the velocity distribution from Delsemme and 
Miller (1971) expressed in terms of a/a max . Using the evaporation rate as re- 
ported by Marsden, Sekanina and Yeomans (1973) , at 0. 8 AU Z — 4. 9 x 10 1 7 and 
a max ~ 15. 5 cm. Since Z varies roughly as r -2 , larger particles can be ejected 
at r <0. 8 AU. Also, as can be seen in Figure 3, at smaller heliocentric dis- 
tances all particles tend to have nearly the same terminal velocity, while at 
larger distances v ~a~^ 2 (approximately). If, as is likely. Comet Encke is 
not covered with ices, then the effective Z(r) should be reduced, thereby also 
reducing a max and v(a). 

For calculating sunlight scattering from dust, we divide the particle size dis- 
tribution into six ranges and assign appropriate mean radii to each range. The 
root-mean-square, (<a 2 > ) 1/2 , is used to calculate the scattering by each size 
range. While (<a 1/2 >) 2 is used for determining the terminal velocity, the re- 
sult is not sensitive to the particular kind of averaging which is used because the 
velocity is a slowly- varying function of particle radius. The size ranges and 
derived properties are given in Table HI-2, with dM/dt = 5x 10 4 g sec -1 and 
a max = 15.5 cm as is appropriate for a heliocentric distance of 0.8AU. 
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Figure HI-3. Terminal Velocities of Non-Volatile Particles 

The first five size ranges used in the Mie scattering code are discussed below; 
the final group contains particles sufficiently large that the geometrical approx- 
imation is adequate ("debris 1 *). The scattering contribution (column 4) repre- 
sents the relative total geometric cross section of each size range. 

Because the particles attain their terminal velocities very near the cometary 
nucleus, throughout the coma (where they will be observed) their spatial distri- 
bution is given by 

P i 

n=(R) = i— cm " 3 

1 4jtR 2 vj 

where Pj is the production rate and v 4 the terminal velocity for particles in the 
i th size range; R is the radial distance from the nucleus. 

This model shares with all current theoretical models the lack of asymmetry in 
ejection which must be present in order to produce the observed non-gravitational 
forces. The size distribution of dust particles lias the virtue of being related to 
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Table E3-2 

Particle Sizes and Terminal Velocities (0.8AU) 


Range 

(pm) 

Production Rate 
(sec -1 ) 

(<a 2 >) 1/2 

(pm) 

Geometric 
Scattering 
Contribution (%) 

V 

(km sec -1 ) 

. 

1 - 1.5 

1.51 x 10 14 

1.28 

9.8 

0.53 

1.5 - 3 

1.79 x 10 14 

2.04 

29.4 

0.49 

3 - 6.7 

3.68 x 10 13 

4.10 

24.5 

0.41 

6.7 - 14.4 

5.03 x 10 12 

9.35 

17.4 

0.33 

14.4-44.4 

8.77 x 10 11 

20.9 

15.2 

0. 25 

44 * 4-a max 

2.34 x 10 10 

62.9 

3.65 

0.15 


observations of actual comets, although there has been no analysis of an old 
comet (such as Encke) using the procedures described by Finson and Probstein 
(1968). Perhaps the lower curve in Figure HI-1, giving Whipple's (1968) dis- 
tribution of interplanetary particles, would be a better representation for Comet 
Encke or the compilation by Dohnanyi (1972). The distribution used here con- 
tains a larger fraction of small particles (1 < a < 5 pm) which are relatively 
efficient scatterers, so our results may represent an upper limit to the sunlight 
scattering for a given dM/dt. 


IV. TRAJECTORY 

Comet Encke is a periodic comet with a period of 3.3 years. It shows greatly 
reduced activity after perihelion relative to before, and shows very little con- 
tinuum radiation at any time, indicating a low dust content (Taylor et al. , 1973). 
Encke never shows a type II (dust) tail and in some apparitions has shown no tail 
at all. The observations from 1885-1951 indicates that Encke's coma becomes 
observable at about 1. 5 AU at the same time as the first appearance of its tail 
(Vsekhsvyatskii, 1964). The orbital characteristics of comet Encke are given 
in Table m-3. 

Vsekhsvyatskii (1964) has estimated a nominal value of 10 5 km for the diameter 
of Encke's coma reduced to a Sun-comet distance of 1 AU. Observational data on 
tail length are given by Yoemans (1973) and TRW (1972). 
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Table m-3 

Orbital Elements of Comet Encke 


Orbital period (year) 

3.30 

Aphelion distance (AU) 

4.09 

Perihelion distance (AU) 

0.34 

Orbital inclination (deg) 

12.4 

Velocity at 1 AU (km s _1 ) 

37.1 

Velocity at perihelion (km s -1 ) 

69.9 

Orbital eccentricity 

0.847 


Many different lands of missions to Encke have been proposed including fast and 
slow flybys. They do have certain elements in common. A nominal miss dis- 
tance is about 10 2 -10 3 km. One proposed probe trajectory calls for a rendevous 
at a heliocentric distance of about R s =0.8 AU. We have chosen this distance to 
calculate the absolute and relative luminosities of the dust, debris, and nucleus. 


V. MIE SCATTERING 

In this section the computer code used to calculate the Mie scattering of sunlight 
from the cometary dust is described. Results are then presented for scattering 
of dust particles in the coma of Comet Encke. 

A. Mie Theory and Computer Code 

Particulate radiation effects in particulate clouds can be calculated provided that 
the optical constant of the radiating particles are known. Gal and Kirch (1973 
and Gal, 1974) reported a new computer code called GMIE to calculate scattering 
cross sections of interest. Calculations are based on MIE theory for a specified 
particulate column with a given complex index for refraction, range of wave- 
length, temperature, and particle sizes. 

It should be noted that the GMIE code can be used to obtain scattering cross sec- 
tions or differential cross section for any given spherical particle or particulate 
column with a given size distribution and index of refraction. 
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1. Electromagnetic Radiation Scattered by Absorbing Sphere 

The passage of electromagnetic radiation through a particle cloud is generally 
accompanied by removal of a fraction of the energy from the incident beam. 

This fraction may be partly absorbed within the particles and partly scattered 
~ i. e. , reappear in the same direction as well as in other directions. The 
characteristics of the scattered radiation are determined by the wavelength X of 
the incident radiation, the refractive index (m = n j - in 2 ) of the particles and 
size as well as the shape of the discrete particles in the medium. For the radi- 
ation studies, we assume spherical particles with radius r. 


The interaction of an electromagnetic wave with an absorbing sphere is described 
by the MIE theory (Fig. III-4) and is discussed in detail elsewhere (van de Hulst, 
1957). Classical MIE theory gives the intensity I(W/cm 2 ) at a distance R and 
angle 9 of the radiation scattered from a single spherical particle of radius r 
exposed to parallel monochromatic radiation of intensity I 0 : 


1 = S(0) 

I 0 2(2?r/X) 2 R 2 


(1II-D 


where S ( 9 ) is the scattering function per unit particle. 


A computational scheme and computer code were provided by J. Dave of IBM 
(1968). These are based on the Ricatti-Bessel functions and give the MIE scat- 
tering cross section of absorbing spherical particles. 


Computations are valid for all values of size parameters that occur in the theory 
regardless of whether they are large or small. The index of refraction of the 
particles are given in terms of their real and imaginary parts. 



Figure III-4. Definition of Single-Particle Scattering 
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The extinction and scattering cross section are calculated from various combina- 
tions of the sum and products of the coefficients a n and b n . The usual expres- 
sions for a n and b n are: 


a 


n 


b 


n 


^n(y) ^ n (x)-m ^ n (y) ^(x) 

(111-2) 

V n (y) £ n (x)-m^ n (y) £[j(x) 


m ifi' n (y) 4/ n (x) - 0 n (y) ^„(x) 

(111-3) 

m i^(y) £ n (x) - ^ n (y) £^(x) 



where 


ra 


= the complex index of refraction 
= a positive integer 
& and £ = the Ricatti-Bessel functions defined by 


n 


(\ v /2 

^ n ( z > = * z j J n + 1 

£(z) = (~n zj Htfl l 


!2 


/ 2 <*> 

(z) 


(IH-4) 


(III-5) 


with J n + l/2 and x/2 the Bessel functions of first and third types. The prime 
denotes differentiation with respect to the argument of the function and 


x = 2ntj\ 
y = mx 


(IU-6) 

(HI-7) 


where r is the particle radius and X is the wavelength. 

Once these scattering coefficients are defined, then extinction and scattering 
normalized cross section can be calculated with the following equations: 


Qext 



(2n + 1) Real (a n + b n ) 


2 °° 

Qscatt (2° + 1) (i a n l 2 + 

x n = l 


(III-8) 
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Scattering efficiency (i.e. , normalized scattering cross section) given by Equa- 
tion (III-9) can also be written as a functi on of scattering angle and the intensity in 
that direction: 


1 r n 

Qscatt =— I [i,(0) + i 2 (0>] sin0d0 (IH-10) 

^ Jr \ 


where 

i,(0) = |S, (0)i 2 

i 2 (d) - is 2 (f?)j 2 


and 


i,(0) + i 2 (0) = S (0) 




ij (0) and i 2 (0) refer respectively, to the intensity of light polarized perpendicu- 
larly and parallel to the plane through the direction of propagation of the incident 
and scattered beam. These intensities are given in terms of the complex ampli- 
tude functions Sj<0) and S 2 (0): 


and 


S,(0) 


v - '' 2n + 1 


[a n 7T n (cos 0) + b n r n (cos 0)] 


(III-12) 


S 2 (0) 


z 


n = 1 


2n + 1 
n(n + 1 ) 


[b n 7 r n ( c ° s 0)+ a n r n (cos 0)] 


(III-13) 


The phase functions 7r n and r n appearing in Equations (III-12) and (III-13) can be ex- 
pressed in terms of the Legendre polynomials, P n , as follows: 


7T n (C0S 0) 


uP n (C0s o) 
d cos 0 


d 7 r n (cos 0) 

r n (cos0) = cos 0 7T n (cos 0) - sin 2 0 

n n d cos 0 


(111-14) 


The detailed computational method to obtain these Legendre polynomials with 
recurrence relationships is discussed elsewhere (Dave, 1968). 
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Otherwise, 


If there is no absorption, i.e., n 2 = 0, then Q EXT = Qscatt • 

Q abs , the efficiency factor for absorption, is given by 

Qabs ” Qext “Qscatt (III-15) 

Equation (EII-9) or (HI-10) gives the total scattered intensity. On the other hand, one 
may need to know the fraction of energy scattered into the forward and backward 
directions. This is obtained in terms of normalized cross section given in Equa- 
tions <III-16) and (ni-17), respectively 

i r^ 2 

Qtrans = ~ T I [i!(0) + i 2 (0)] Sin0 dO (111-16) 

* J o 

and 

l r n 

Qref = 7 CijC^) + 1 2 (0)] sin0d0 (111-17) 

X Jl 7/2 

2. Solar Scattering Cross Section for a Spherical Particle 

With the help of the basic MIE scattering theory, particular scattering due to 
sunshine will be treated through an equivalent normalized differential solar- 
scattering cross section. The sun is approximately half a degree as observed 
from the earth, and calculations similar to those derived by Gal and Kirch (1973) 
for earthshine scattering would require evaluation of scattering function S(0) for 
smaller than 1 deg, which is our present minimum step size. Approximate dif- 
ferential cross section will be obtained by averaging the basic MIE scattering 
function over a small angular increment (A 8 ) to account for the finite solar disk 
(Fig. m-5). 



Figure m-5. Solar Scattering Model 
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From Equation (131-1) 


(I S R 2 ) 1 



m 


(III- 18) 


The normalized differential cross section can be obtained by dividing Equation (in-18) 
with the geometrical cross section G = r 2 ir 

( I R 2 \ 1 

_L_] = S(0) (sr 1 ) (111-19) 

" 2 .( 3 =) 

Equation (m- 19) yields the normalized scattering cross section for a single spherical 
particle. 

3. Column Average Cross Section and Radiance 


a. Column Scattering — The extension to the case of scattering from a column 
of particles is straight-forward upon adopting restrictive assumptions. First, 
the distance R must be large compared to the column dimensions so that scatter- 
ing angle is essentially constant for a given I 0 direction. Also, the column must 
be optically thin. This requirement means that: (1) each particle is exposed to 
the incident intensity I 0 , regardless of location; and (2) the radiation scattered 
from one particle does not interact with others as it passes from the column; 

i. e. , there is no multiple scattering. 

These assmnptions allow the scattered intensities to be added. With the help of 
Equation (in-19) we can define a cloud averaged differential cross section, a. 

a = X ( d Q)j A *(Nj/N t ) (pm 2 sr 1 ) (111-20) 

i 

where 

N t - the total number of particles 
Nj = the number i th particles 
r . - the radius of the i th particle 

b. Column Radiance — Figure H3-6 shows the geometry; the solid angle of the 
instrument is given by £2 = A/L, and the scattering cloud volume is V = AC, 
where A is the column surface area, and £ is the depth of the column. The total 
number of particles in the scattering values is N t = nAC, where n is the particle 
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Figure in-6. Definition of Scattering Values 


density. The incident intensity I 0 x (watt/cra 2 • p) can be calculated assuming 
the sun as a blackbody at a temperature T sun = 6000 K, and sun is approximately 
half a degree as observed from the earth 

I oX = B X (T = 6000) (Aa) (Won" 2 pm’ 1 ) (111-21) 

where 

1.19 10 4 (1/A) 5 

— (W cm 2 sr 1 pm l ) (111-22) 

( 1.4384 X 10 4 \ 

where A is the wavelength in microns, and (Ax) is solid angle of the sun viewed 
from the column (6.8 10 -5 sr) at 1AU. 

The column radiance is defined as 

N x = — (W cm -2 sr -1 pm _i ) (III-23) 

A 

Substituting Equations (III-19) through (III-22) into Equation (III-23), we obtain 

5 

Nx = Bxn T l(Ao:) ^ (dQ)j r 2 7r (Nj/Nj) (W cm -2 sr -1 pm -1 ) (111-24) 
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Equation (III- 24) has been programmed into our GMIE code; additional inputs re- 
quired are particle size, its distribution, and the column depth, £. 

B. Mie Scattering From Cometary Dust Particles 

We have inputed the comet Encke dust distribution into the GMIE code for vari- 
ous indices of refraction. The normalized scattering and absorption cross sec- 
tions given inEquations (III-10) and (HI-15) are plotted in Figure HI-7 for an index of 
refraction of 1. 7-0. 05 i. Figure III-8 shows a plot of the scattering function as 
a function of scattering angle for a dielectric particle of 0. 65 pm radius with an 
index of refraction n = 1. 7-0. 05 i. More plots of particle scattering functions 
are found in Appendix I showing the effect of particle sizes and indices of refrac- 
tion. Also, the angular dependence of the polarization is plotted for a few select 
cases. To calculate the scattered sunlight from the comet, we chose an index of 
refraction of 1.7-0. 05 i to represent cometary dust particles which are probably 
primarily silicates containing metallic elements (Wiekramasinghe and Krishna 
Swamy, 1968), Column densities are calculated in Appendix H. The results of 
the Mie scattering calculations are presented in Table DI-4 where the wavelength 
dependence of the scattered radiance N x from Equation (III-23) is presented as a 
function of the Sun-Comet-Probe angle. The functional dependence on wavelength 
and angle of the average differential cross section a from Equation (IH-20) is dis- 
played graphically in Figures in -9 and EEI-IO for a few selected cases. Figure 
HI-9 shows the effect of view angle on a at A =0. 5pm. The minimum scattering 
at this wavelength can be seen to occur in the 0-90° region. The wavelength de- 
pendence of the scattered radiance is displayed in Figure HI-10 for several view 
angles. It would appear that the scattering in the 0-90° range is relatively in- 
sensitive to wavelength dependence. 

These results, which have been calculated for a heliocentric distance of 
R s -0.8 AU, can be approximately scaled to other distances by the law 

( 0.8\ 4_ f 

N X (R S ) = N x (0.8) J (IH-25) 

where f accounts for the terminal velocity dependence on heliocentric distance 
and R s is inAU. Equation (HI-25) reflects the inverse square dependencies of the 
solar irradiance and dust production rates. The heliocentric dependence of the 
terminal velocities arises primarily in the case of Encke from the fact that the 
acceleration of the dust particles is dependent on the number of collisions the 
particles undergo with effusing gas molecules. Hence, the terminal velocities 
depend on particle sizes and the vaporization rate of gas from the nucleus. In 
the case of Encke dependence on the temperature of the nucleus is of lesser 
importance. For the particle sizes of importance here the heliocentric depend- 
ence can be approximated by setting f - 0.4. This law will hold for R s out to 
about 1.5 AU where the dust production ceases due to the threshold effect of ice 
evaporation. 
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Figure III-8. Angular Dependence of Scattering Functions, a = 1.28/xm 

n = 1. 7-0.05 i 
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Table HI-4 


Radiance of Comet Dust at R s =0.8 AU 
N^(w cm -2 sr -1 ^nf 1 ) 



Sun-Comet- Probe Angle 

0° 

30° 

60° 

90° 

120° 

150° 

180° 

0.3 

2. 7(-6) 

2. 4 (-6) 

2. 3 (-6) 

2, 2(-6) 

6. 2 (-6) 

1. 2 (-5) 

1.4 

0.4 

4. 8 (-6) 

4. 2 (-6) 

4.1 (-6) 

4.1(-6) 

8. 3(-6) 

2. 3(-5) 

1.5 

X0.5 

5. 4(-6) 

4. 6(-6) 

4.4(-6) 

5. 4(-6) 

8. 9(-6) 

2. 6 (-5) 

1.1 

(fiin) 0 . 6 

5. 0(-6) 

4 . 1(— 6) 

4.1(-6) 

6.0 (-6) 

9.0 (-6) 

4. 2 (-5) 

6. 8(-l) 

0.7 

4. 3(-6) 

3. 5(-6) 

3. 3 (-6) 

4. 7 (-6) 

6. 8 (-6) 

9. 9(-6) 

4. 2(-l) 

0.8 

3. 7 (-6) 

2. 9 (-6) 

2. 5 (-6) 

3. 6(-6) 

6. 3 (-6) 

2. 8 (-5) 

2. 7 (-1) 

0.9 

2. 9(-6) 

2. 3(-6) 

2. 2 (-6) 

3. 4(-6) 

4. 0(-6) 

3.0 (-5) 

1.8(-1) 

1.0 

2. 6(-6) 

1.8 (-6) 

1. 8 (-6) 

3. l(-6) 

1.4(-6) 

1.8 (-5) 

1.2(-1) 


VI. SCATTERING BY DEBRIS PARTICLES 

For particles in the 10" 2 -10 cm range, scattering can be treated essentially as 
isotropic with a geometrical cross section. Thus, the differential cross section 
can be written as 


da 7ra 2 a n 
dfl 47r 


(111-26) 


where the albedo a n is taken equal to that of the nucleus, 0.1. Calculating an 
average differential cross section for R s = 0. 8 AU using the distribution in 
Table HI-2 results in 


da a 2 a_ 

= — — = 9.9 X 10 1 fxm 2 sr -1 (111-27) 

dO 4 
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where a 2 is the rms radius of large particles, a > 44 pm. The column density 
for the debris particles is given by Equation (B-3) which results in 

P= 1 

n T = = 2.5 X 1CT 1 particles cm 2 (111-28) 

v i 47rR min 

Using Equations (in-27), (III-28) with Equation (HI-24), we then would get a scattered 
brightness for large particles at R s = 0. 8 AU for X = 0. 5 pm of 

N x = 8.9 X 10~ 8 W cm" 2 sr" 1 pm" 1 (111-29) 

As can be seen from Table HI-4, this contribution is less than the small particle 
Mie scattering by a factor of 50. So it can be concluded, following our model, 
that the small size dust particles contribute much more to the scattered bright- 
ness than the large size debris particles. 


VH. SCATTERING BY ICE PARTICLES 

It is difficult to precisely assess the ice particle distribution and production rates 
for Encke. We shall assume here that the ice particles can be represented by 
particles with a radius of 300pm (Delsemme and Miller, 1971). An upper limit 
for the ice production rate can be derived from the OGO-5 satellite observations 
of Lyman-a emission from Encke. The observed emission corresponded to a 
production rate of Q(H) = 5 x 10 26 atoms sr -1 s -1 for hydrogen atoms at R $ = 
0.715 AU (Bertaux et al. , 1973). Assuming that all the H atoms came from ice 
particles results in a production rate for water Q(H 2 0) = 3 x 10 2 7 molecules s -1 
(Delsemme and Rud, 1973). For icy grains of 300 pm radius, we then have 
Q ice ~ 2 x 10 9 s' 1 . 

Using these numbers, it is possible to estimate the scattered sunlight assuming 
an albedo of 0. 9. For the differential cross section we have 

= 2.0 X 10 4 pm 2 sr" i (III-30) 

dil 

and the column density would be 

n T ~9 X 10 -2 particles cm -2 (III-31) 

An estimate of the scattered light would then yield for X = 0. 5 pm and R s =0.8 AU 

N x = 6 X 10~ 5 W cm" 2 sr" 1 pm' 1 (111-32) 
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This result is of the same order of magnitude as the results for Mie scattering 
for cometary dust suggesting that scattering by ice particles surrounding the 
nucleus could be a major contributer to scattered sunlight. On the otherhand, 
our ice model is tentative and ice production rates could be revised downward. 


Vin. SUNLIGHT REFLECTION FROM THE NUCLEUS 

We have assumed an albedo of 0. 1 for the nucleus so the brightness of the nucleus 
following Equation (HI-24) is given by 


N£ uc = B x (Aa) (III-33) 

where Aa is the solid angle subtended by the sun at a given heliocentric distance, 
R s . It varies as R~ 2 . It is assumed here that the scattering is isotropic into 
2irsr. For X = 0.5 pm and R s =0.8AU, we have 

N£ uc = 5.7 X KT 3 WcnT 2 sr 1 pm" 1 (III-34) 

IX. COMPARISON OF BRIGHTNESS OF NUCLEUS TO BACKGROUND 
BRIGHTNESS 

The contribution to the background brightness from the comet dust is- tabulated 
in Table HI-5. Taking one specific case we can compare the brightnesses from 
the nucleus to the background brightness. For the case with a heliocentric dis- 
tance of 0.8AU, A = 0.5 pm, and view angle a = 90°, the radiance levels are 
presented in Table III-IV. 

The background from particulate scattering is 2 to 3 orders of magnitude below 
the intensity from the nucleus. The relative particle to nuclear brightness levels 
in Table m-4 are expected to vary less than an order of magnitude going into 
perhelion (0. 34 AU) or out to "turn on" at about 2 AU. 


Table HI-5 

Radiances at R s = 0.8 AU for A = 0. 5 pm, a = 90 9 


Nucleus 

Particulate Scattering 

Dust 

Debris 

Ice 

Radiances 5.7 x 10 -3 
(W cm -2 sr -1 pm -1 ) 

5.4 x lO" 6 

8. 9 x 10“ 8 

6 x 10' 6 
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Based on these calculations there should be no difficulty in viewing the nucleus from 
masking by particulate scattering. The particles are optically thin and the bright- 
ness levels are several orders of magnitude down from that of the nucleus. 

X. OPTIMIZATION OF IMAGING SYSTEM 

In terms of scientific value, resolution as high as possible is desired to examine 
the nature and composition of the surface. One of the problems with high resolu- 
tion imagery and perhaps the limiting one is the large relative velocity of the 
spacecraft and comet which can be of the order of lOkm/sec. This, of course, 
would depend on the particular mission trajectory selected. Blurring of the 
image while a pixel cell is being exposed necessitates short exposure times. 

This requirement would argue for a high sensitivity-low noise system such as an 
image disector or a spin-scan camera. As an example of some numbers con- 
sider the case where the desired resolution element is ten meters squared on 
the surface of the nucleus and the bandwidth is 0. 1 ji. Then the photon rate cor- 
responding to a resolution element would be N = 1 x 10 21 photons-sec" 1 -sr -1 . 

The solid angle subtended by a 5” diameter telescope at 1000 km would be about 
2. 5 x 10 -1 5 ster and would result in photon rate at the detector of 3 x 10 6 photons/ 
sec per pixel cell. If the transverse velocity of the comet with respect to the 
spacecraft is of the order of lOkm/sec then exposure times/pixel of no more 
than 1msec are required. This corresponds to 6 x 10 3 photons/pixel which is 
sufficient to detect intensity variations from the surface of 4%. It should be 
noted that the intrinsic readout noise of a vidicon system corresponds to about 
1000 electrons/pixel. 

Other options which should be considered with an imaging system are use of 
polarizers to gain polarization information about the nuclear surface. Because 
the illumination of the nucleus will be quasi-lambertian, viewing at different 
phase angles will afford information regarding surface irregularities. Shadow 
heights can be used to topographically map the surface. The use of filters to 
block out line and band emission features in the bandwidth of the imaging system 
may be deemed wise. 

It is concluded that dust and debris scattering will not compromise the nucleus 
imaging mission and that high relative velocities warrant the use of a low noise 
imaging system such as a spin-scan camera. 
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Figure A5. Angular Dependence of Polarization, a = 0, 65 pm, n = 1. 5-0. li 
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CALCULATION OF COLUMN DENSITIES 


The number of particles N passing through a spherical surface at a distance R 
from the center of the nucleus during a time At is given by 

N = p i At = £ nj(R) vj At 4 ttR 2 (B-l) 

i i 

where P i5 nj and Vj are the production rate, the number density and the terminal 
velocity, respectively, of the particles in the ith size range. 


The number density as a function of R is then given by 


nj(R) = 


P i 

Vj 4irR 2 


(B-2) 


The column density of particles n co ] is obtained by integrating Equation (B-2) 


from R mln to R 


min 


max 


v-i r R ma 
n col = 2^ / 

j jR min 


max PjdR 


4tt VjR 2 


(B-3) 


R m in is the radius of the nucleus (R n ), and since Vj is constant for R > 10 R,,, 
the value of n co] is insensitive to the value of R max . Equation (B-3) is actually 
the lower limit to n co j , since for R < 10 R n the particle velocities are less than 
Vj. The calculations of Delsemme and Miller (1971), shown in their Figure 2, 
suggest that the particles are subject to approximately constant acceleration un- 
til they move at nearly terminal velocity. The integral in Equation (B-3) is most 
sensitive to the particle velocities at small R. As an upper limit to n col , sup- 
pose the actual particle velocities are given by 


v(R) 


v o R n <R<10R n 


^ R> 10 R n 


(B-4) 


Here v 0 is the initial velocity at which particles are ejected from the nucleus. 
The column density for species i is then 

Pj 

"col = [0.1 +0.9(Vj/v o )] (B-5) 

4tt R nVi 

Taking v Q = 0. 1 Vj gives a column density 9. 1 times the value obtained from 
Equation (B-3) assuming v = Vj everywhere. We could, of course, go on to 
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evaluate Equation (B-3) using the assumption of constant acceleration to obtain 
somewhat better estimates for n col , but the results depend on the unknown ejec- 
tion velocity v Q . For the present, we simply point out this refinement which 
should be included in more comprehensive models. 

For calculating impacts suffered by a spacecraft passing near the comet, we 
have evaluated the column density for a linear path passing within R Q km of the 
nucleus (we will call this the impact parameter). The complications discussed 
above are not encountered because R 0 5t> R n . The spatial density is given by 
Equation (B-2) , and the column density is 


a col^ 4 ’ R o) 


P(a) 

4v(a)R 0 


(B-6) 


Since v(a) -*> 0 as a a max , the production rate P(a) must also go to zero in 
order to avoid a singularity. Since the terminal velocity for a > 10 “ 4 a max is 
given by 


v(a) - 2.88 X 10 3 (a/a max ) 1/2 [1 - (a/a max ) J 12 ] km/sec (B-7) 

(from Delsemme and Miller, 1971) , we cut off the production rate by a factor 
[1 -(a/a max ) 1,r2 ] . This does not affect the results presented in Section III. The 
production rate for a > 4.44 x 1(T 3 cm is given by 

P(a) = 7.298 X 10~ 4 M dust a' 5 [1 - (a/a max ) 1/2 ] sec" 1 (B-8) 

where Mdust — 5 x 10 4 g sec -1 and a max =15.5 cm at heliocentric distance 
0.8AU, while Md ust — 2.77 x 10 5 g sec -1 and a max = 98.2cm at 0.34AU 
(perihelion) . 

Using Equations (B-7) and (B-8) in Equation (B-6), then integrating over particle 
sizes from a to a max , we obtain the cummulative column density: 


n(a, R 0 ) 



(B-9) 


M 

= 1.813 X 10- 7 - - ^ ust - [a" 3 5 - a ^ 3 ax ] (B-10) 

a max R o 

this represents the number of particles which have radii a or larger. Table B-l 
gives representative values. 
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CUMULATIVE COLUMN DENSITY (CM 



PARTICLE RADIUS (CM) 

Figure B-l. Cumulative Column Density vs. Particle Radius 
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Table B-l 

Column Densities For an Impact Parameter of 100 km 


Heliocentric 

n{a, 100km) (cm -2 ) 

Distance (AU) 

a = 10" 2 cm 

a = 10 1 cm 

a = 1 cm 

0.34 

5.1 x 10- 3 


1.6 x 10“ 6 

5.1 x 10" 10 

00 

O 

2.3 x 10" 3 


7.3 x 10" 7 

2.3 x 10" 10 


For given values of a and R Q> n varies approximately inversely as the heliocen- 
tric distance because M dust 311(1 a max ^ otl1 vaI T approximately as r~ 2 . 
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APPENDIX IV 


MODEL OF IMAGING EXPERIMENT FOR 1980 ENCKE MISSIONS 


A complete description of the photometric characteristics of comet Encke in- 
cludes three targets : coma, halo, and nucleus. 

A. Coma 

The model representing Encke's coma appearance at 1 AU comet-sun distance 
has been given by F. Taylor (Ref. 2). A relative brightness directly propor- 
tional to angular distance from the coma center has been combined with visual 
observations of the central brightness (B Q ) to yield a function applicable to 
within 10 5 km radius of the nucleus. At comet-sun distances (A) other than 1 AU, 
the relative brightness is scaled according to the empirical observation 

B 0 = Bo x lO -0 - 4 * 15 Io g A) A > 0.7 AU 

The signal recorded by a camera of given focal length is a function of both 
pointing within the coma and area of the projected pixel field -of- view. For this 
reason, it si convenient to compute the exposure obtained at each element of the 
camera's sensing device (Ref. 3): 


where 


1 7T f d / 2 r n 

Exp. - E 0 / dw — — I — 27rr dr 
4T Jo r + r o 


E 0 / dw “ illuminance/str at the coma center 

t / 4 T 2 = camera field-of- view with optics transmission and shutter speed 

d = FOV pixel x to comet) 

/ dr = area of projected pixel convolved with brightness function 

Because the intensity of coma emission is not radially symmetric over a given 
pixel, the following approximation is used: 

Exp = BMt + 

JTX 2 


rv-i 



where 


E(x) 


E 0 /dw 

* 7 r 

4T 2 


= illuminance at pixel n, distance from coma center (x n = nd ±d/ 4) 
Ax = x 2 - x 1 

Finally, the light from the coma consists of discrete spectral emission bands. 
The total integrated brightness within the sensitivity of the camera is given by: 


= E E n S X T * 
n 


where 


n = spectral line 

= camera sensor spectral sensitivity 
T x = telescope spectral transmission 


For example, according to the Taylor model, the coma illuminance per steradian 
at 1. 2km from the nucleus is 0. 36 fc/str due to the C 2 line at X = 5165 A. In- 
tegration of all spectral lines produces an effective energy for GEC vidicons of 
1.4 fc/str. The corresponding exposure at pixel n obtained by an approaching 
spacecraft depends only on the optics-sensor configuration (4T 2 , FOV pixel), 
the spacecraft-comet range, and the comet-sun distance. 


10 (5 - 6 log A) (Ax+ i.i) (A log (x+ 1.1)) 
Exp p = 7 jrt 


4T“ 


x 2 


where 

d = FOV pixel x Range (km) 

A = comet-sun distance (A U) >0.7 AU 
t = shutter speed (sec) <1 pixel S 
x = nd-d/4 
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X 


= ( Xl + x 2 )/2 
t = 1.3 t/2 2 pixels smear 
t = 1.9 t/5 5 pixels smear 
e.g. , Spr/sec 1 pixel smear = 6 sec 
B. Halo 


As discussed in the Taylor report, a continuous background spectrum is also 
observed in the spectra of Encke. At 1/5 the E Q /dw strength of the C 2 (X5165A) 
emission light, this reflected sunlight suggests the presence of a halo to a rad- 
ius of approximately 1000 km. Regardless of whether this cloud is composed of 
icy grains or dust, an aerosol-scattering model (conservative) may be fitted to 
those observations. At comet-sun distances other than 1AU, both the reflected 
light and the halo optical depth should be scaled proportional to the square of the 
solar distance: 


where 


E oh /dw = (E Qh /dw) 1AU x 


1 - e 


,-tM 


L 1 - e - toM 


) x ^ 


x wo 


t = optical depth (t Q A~ 2 ) 
to = depth at 1 AU ~0.1 (E Q /dw = 1/5 C 2 at 1 AD) 
wo = scattering indicating ~1. 0 


M = air mass 


- 2ati = e = g = 0 


Again, fitting an r _1 brightness function to r = 10 3 km (compromise between icy 
grains and dust particle scattering) the exposure is calculated as follows: 


Exp 


ffl5t(x + r 0 ) (log (x + r Q ) 
4T 2 x 2 




r < 1 0 3 km 


C . Nucleus 


A variety of possible nucleus sizes and albedos are predicted by differing bright- 
ness observations (which gives the product of albedo x area of surface disk). 
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Roemer's (Ref. 2) photography of Eneke at aphelion (minimal coma activity) sug- 
gests pr 2 = 0. 8; e.g. , r = 2km, p = 0.2. Other observations over a large frac- 
tion of the orbital path indicate pr 2 >0.24 (r = 2km, p = 0.06) or an apparent 
magnitude: 


Mj * 16.0+ 5 log A+ 5 log R + 0.03 g 


where 

R = range to comet (AU) 

0.03 g = lunar phase function 
g - phase angle (deg) 

Calculations in this memo use only M 7 . The difference in intrinsic brightness, 
for example, at R = 2. 3 x 10 6 km, A = 0. 8 AU, g = 13° is: 

pr 2 < p 0 1/2 

M. = -26.4 -2.5 log = 26.15 -5 log —— - 5.6; pR 2 = 0.8 

1 a 2 r 2 a r 

M 2 = 6.8, pR 2 = 0.24 

where 

pr 2 = 0.8 
R = range (km) 

A = (AU) 


The camera response to a point source is predictable from television readout 
theory or by reference to flight experience: 

Exp = 10- (6 - 71 + °- 4M) t — TFf 

a 

where 


10 -6 - 71 = illuminance (fcs) from a 0 th mag star 


— = ratio of aperture to image area 
a 

T = optics trans m ission 
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F = filter transmission 

f = television read beam efficiency factor, f may be determined by 
Mariner 9 star photography 

For example, the Mariner 9 B-camera with clear filter and MVM optics IR coat- 
ings (T = 0.525) receives an exposure of: 

Exp = 10.8 X 10-°- 4M t 

As an extended source, the sensor exposure is given by: 


Exp 


M e -.M Jl 

A 2 4T 2 


where 


S = solar constant ~13,000fc 


The above camera would yield a signal-to -noise response (S/N = 10 at O.Olfcs) 
to the nucleus of Encke (p = 0. 2, g = 0° , A =0.54 AU) : 


9000 

S/N = 10 X t F e -tM 

4T 2 


0.01 fcs 


10 3 X 110 tF 


The ways on which the expected imaging results are affected by the intercept 
mode are shown in Tables IV-1 and IV-2 and in Figures IV-1 thru IV-8. 
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Table IV-1 

Comet-Sun Encounter Range Effects on Imaging Experiment* 


Parameter 

Encounter - 
Sun Distance 


0. 34 AU 

0.53 AU 

0. 80 AU 

1. Flyby velocity at encounter (km/sec) 

8.6 

18 

26.8 

2. Comet Observation Potentials 




Coma maximum signal to noise ratio at 
E -10 d (clear filter, 6 second exposure) 

45 

18 

6 

Earliest halo detection (Range - km) 

E -14 d (10 7 ' 1 ) 

E -3 d (10 6 - 7 ) 

E -l d (10 6 - 4 ) 

Earliest nucleus resolution (Range - km) 




M = 16 + 5 log R + 5 log + 0. 03 g: 
Without coma (10% filter): 

E -3 d (10 6 ) 

E -7 d (10 6 ' 8 ) 

E -1. 6 d (10 6 - 4 ) 
E -2.5 d (10 6 - 6 ) 

E -0. 8 d (10 6 - 3 ) 
E -1. 6 d (10 6 - 5 ) 

Smear rate at encounter (pixels/sec) 

287 

600 

895 

Image smear for S/N = 50 with p = 0.2(0.05) 
(pixels) Enc. = 1500 
1000 
500 

0.1 (038) 
0.14 (0.57 
0.28 (1.14) 

0.4 (1.6) 
0.6 (2.4) 
1.2 (4.8) 

1 (3.8) 

1.4 (5.6) 
2.8 (11.2) 

3. Number of possible pictures 




Encounter duration with 2 km > 5 pixels 
(Seconds): 

3000 

1480 

545 

Number of pictures without a tape recorder 
Large frames** (84sec/fr): 

33 (50Kbs) 
4 ( 4 Kbs) 

18 (120 Kbs) 

7 (120 Kbs) 
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Table IV-1 (Continued) 




Encounter - 


Parameter 


Sun Distance 



0. 34 AU 

0.53 AU 

0.80AU 

3. Number of possible pictures (Cont. ) 




Large frames (10 sec/fr): 

33 (50Kbs) 
4 ( 4 Kbs) 

39 (120Kbs) 

14 (120 Kbs) 

Small frames*** (15 sec/fr): 

94 (50 Kbs) 

39 (120 Kbs) 

36 (120 Kbs) 


8 ( 4 Kbs) 



Number of pictures with tape recorder 
(10 9 bit capacity) 




Large frames (84 sec/fr): 

36 

18 

7 

Large frames (10 sec/fr): 

240 

148 

55 

Small frames (15 sec/fr): 

200 

148 

55 

Linear rate per 84 sec (km/picture) 

720 

1510 

2250 

4. Spacecraft/Mission Concerns 

0. 34 AU 

0. 54 AU 

0. 78 AU 

Hardware essentials: 

FEPS, Tape recorder 
X band transmission 

FEPS 

— 

Round trip light time at Encounter 

16 min 

9 min 

5. 5 min 

Significant attributes (S/ C) at 

Fly by velocity 

Compromise 

Earth-S/C 

Encounter 



distance 

Significant comet advantages: 

Nucleus luminosity. 
Coma activity 
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Table IV-1 (Continued) 


Parameter 

Encounter - 
Sun Distance 

0. 34 AU 

0. 54 AU 

0. 78 AU 

4. Spacecraft/Mission Concerns (Cont. ) 
Significant mission problems: 

Automatic sequencing 
(No RT data with 
nucleus as an ex- 
tended source) 

Same 

Same + number of 
nucleus pictures 


*10 3 Km Encounter; 50 cm fl, f/2.35 Vidicon Camera; Despun Platform 
**Large Frame: 700 x 800 elements x 8-bit/element = 4. 5 x 10 6 bits 
***Small Frame: 250 x 800 elements x 8-bit/element = 1.6 x 10 6 bits 
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Table IV- 2 


Comet-Spacecraft Encounter Range Effects on Imaging Experiment* 


Parameter 


Closest Encounter** 

A(r = 0) 

B(r = 0) 

A(r = 400) 

B(r = 400) 

1. Encounter duration (Image 

0. 34 AU 

606 sec 

592 sec 

593 sec 

568 sec 

> 24 pixels): (seconds) 

0. 53 AU 

289 

284 

282 

266 


0.80 AU 

194 

193 

190 

179 

2. Number of pictures possible 

0.34 

30 

29 

29 

28 

at 20 sec/frame (Image S 6 

0.53 

14 

14 

14 

13 

25 pixels) 

0.80 

10 

10 

10 

9 

3, Encounter Image Resolution 

0.34 

6m 

29 m 

27 m 

49 m 

at 0° phase (0.1 modulation) 

0.53 

10 

23 

32 

58 

(midscale exposure, 10% 

0.80 

46 

31 

37 

58 

albedo) 






4. Encounter Image Resolution 

0.34 

32m 

33m 

30 m 

50 m 

at 90° phase (0.1 modula- 

0.53 

79, 

54 

55 

66 

tion) (midscale exposure, 

0.30 

260 

240 

230 

153 

10% albedo) 






5. Angular Encounter Rate 

0.34 

5°/ sec 

0. 9° /sec 

1. 0°/sec 

0. 6°/sec 

(deg. /sec) 

0.53 

10.2 

2.5 

1.8 

1.0 


0.80 

23.9 

4.9 

2.9 

1.5 
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Table IV-2 (Continued) 


Parameter 


Closest Encounter* ** 

A(r = 0) 

B(r = 0) 

A(r = 400) 

B(r = 400) 

6. Data rate per image segment 

0.34 

390 KBS 

73 KBS 

79 KBS 

43 KBS 

at encounter required for 

0.53 

380 

93 

66 

37 

10 second frame time 

0.80 

605 

124 

73 

38 

(Kb/sec): (Image diameter 






x 800 samples x 7 bits/ 






sample) 






7. Image size at Encounter 

0.34 

700 pixels 

131 pixels 

141 pixels 

77 pixels 

(2 km nucleus): 

0. 53 

680 

166 

118 

66 


0.80 

1080 

221 

130 

67 


*50 cm fl, f/2.35 vidicon camera, despun platform 

**0. 34 AU: A(r = 0) = 95 km (optical navigation) B(r = 0) = 510 km (w/o optical navigation) 

A(r = 400) = 471km (optical navigation) B(r = 400) = 865km (w/o optical navigation) 
0.53AU: A(r = 0) = 98 km (optical navigation) B(r = 0) = 403 km (w/o optical navigation) 

A(r ~ 400) =566 km (optical navigation) B(r = 400) = 1017 km (w/o optical navigation) 
0. 80 AU: A(r = 0) = 62 km (optical navigation) B(r = 0) = 302 km (w/o optical navigation) 

A(r = 400) = 515 (optical navigation) B(r= 400) = 996 km (w/o optical navigation) 




Figure IV-1. This figure compares the relative brightness of Encke's coma, 
halo and nucleus image center as a function of spacecraft distance for a 0. 34 AU 
comet-sun distance. With a camera threshold of 1CT 2 fcs (S/N = 10), the broad 
bandpass coma image is detected more than 20 days before encounter; however, 
the nucleus does not exceed background coma emission until E - 3 days. 
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Figure IV-2. This figure compares the relative brightness of Encke's coma, 
halo and nucleus image center as a function of spacecraft distance for a 0. 53 AU 
comet-sun distance. With a camera threshold of 10-2 fcs {S/N = 10), the broad 
bandpass coma image is detected more than 15 days before encounter; however, 
the nucleus does not exceed background coma emission until E - 2 days. 





Exposure q-fc Vud<e*H 



Figure IV-3. This figure compares the relative brightness of Encke's coma, 
halo and nucleus image center as a function of spacecraft distance for a 0. 80 AU 
comet-sun distance. With a camera threshold of 10" 2 fcs (S/N = 10), the broad 
bandpass coma image is detected more than 5 days before encounter; however, 
the nucleus does not exceed background coma emission until E - 1 day. 
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Figure IV-4. The apparent coma image size (abcissa) and 
signal-to-noise ratio (ordinate) are shown in this graph to be a 
function of shutter speed and spacecraft range. At 0. 8AU, as- 
suming a spacecraft attitude drift rate of 5ur/sec, a 30-second 
shutter speed would yield an image nearly twice as bright as a 
6- second exposure but would exhibit 5 picture elements of 
smear. The camera threshold of 10" 2 fcs, however, is not 
quite attained at E - 10 days using a clear filter. 
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IMAGE DIAMETER, pixel* 


Figure IV-5. The apparent coma image size (abcissa) and signal-to- 
noise ratio (ordinate) are shown in this graph to be a function of shut- 
ter speed and spacecraft range. At 0. 53 AU, assuming a spacecraft 
attitude drift rate of 5ur/sec, a 30-second shutter speed would yield 
an image nearly twice as bright as a 6-second exposure but would ex- 
hibit 5 picture elements of smear. The camera threshold of 10' 2 fcs, 
however, is barely attained at E - 20 days using a clear filter. 


IV-15 




IMAGE DIAMETER, pixels 

50 g o-5 



RANGE, kfr 

Figure IV-6. This figure illustrates the required shutter speed to obtain an 
image of the Encke nucleus at S/N^SO. The image S/N is dependent on the 
nucleus brightness which is a function of albedo (0. 05 vs. 0.2), phase geometry, 
and distance if the image is a point source (>10 5 km range). Following the ap- 
proaching nucleus from right to left, it is apparent that the required shutter 
speeds yield less than one pixel of image smear (clear filter) unless narrow 
bandpass filters are used (Filter Factor = 20). 
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SHUTTER SPEED TO PRODUCE $/N 
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Figure IV-7. This figure illustrates the required shutter speed to obtain an 
image of the Encke nucleus at S/N= 50. The image S/N dependent on the nucleus 
brightness which is a function of albedo (0.05 vs. 0.2), phase geometry, and 
distance if the image is a point source (>10 5 km range). Following the approach- 
ing nucleus from right to left, it is apparent that the required shutter speeds 
yield less than one pixel of image smear (clear filter) unless narrow bandpass 
filters are used (Filter Factor = 20). 
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Figure IV-8. This figure illustrates the required shutter speed to obtain an 
image of the Encke nucleus at S/N=50. The image S/N is dependent on the 
nucleus brightness which is a function of albedo (0. 05 vs. 0.2), phase geometry, 
and distance if the image is a point source (>10 5 km range). Following the ap- 
proaching nucleus from right to left, it is apparent that the required shutter 
speeds yield less than one pixel of image smear (clear filter) unless narrow 
bandpass filters are used (Filter Factor = 20). 


IV-18 




